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ABSTRACT 


This  study  is  an  analysis  of  the  performance 
characteristics  of  a  trench  drain  system  used  for 
foundation  dewatering.  By  the  uses  of  a  finite  element 
analysis  program  the  trench  drain  system  performance  will 
be  modeled. 

Through  the  use  of  dimensional  analysis  techniques 
the  results  of  the  system  modeling  will  be  used  to  prepare 
a  means  of  predicting  the  system  performance  for  a  wide 
range  of  variables. 

As  an  economic  concern  the  trench  drain  system  will 
be  compared  with  the  common  blanket  drain.  This 
comparison  will  provide  information  necessary  for  the 
selection  of  the  method  most  economical  for  the 
application. 

The  trench  drain  characteristics  of  primary  concern 
are  the  total  system  flow  and  the  maximum  free  surface 
height  of  the  water  between  trenches.  The  determination 
of  the  radius  of  drawdown  for  the  system  is  a  vital 
element  of  the  prediction  process  and  is  a  major  portion 
of  the  study. 
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CHAPTER 


TRENCH  DRAIN  SYSTEM 

I .  Introduction: 

Shallow  foundations  are  susceptible  to  damage  and 
leakage  from  the  hydraulic  forces  of  groundwater.  The 
pressure  of  water  will  damage  walls,  pass  into  the 

structure  through  cracks  and  holes  and  contribute  to  the 
deterioration  of  metal,  wood  and  concrete  structural 
members.  A  structure  will  not  satisfactorily  full  fill 
the  intention  of  the  owner  if  water  gathers  in  the 

basement  or  garage  of  the  building.  A  great  deal  of 
effort  is  expended  by  builders,  owners  and  designers 
attempting  to  prevent  water  from  infiltrating  a  foundation 
or  collecting  behind  the  walls.  The  most  simple  solution, 
as  well  as  the  most  effective  and  economical  is  to 

properly  design  the  foundation  so  that  water  is  not 

present  to  penetrate  the  structure. 

For  a  shallow  foundation  positioned  at  a  moderate 
depth  below  a  water  table  the  foundation  trench  drain 
system  is  an  alternative  to  expensive  water  resistent 
methods  of  construction.  The  foundation  trench  drain 
system  is  an  application  of  a  very  old  technique.  A 
trench  filled  with  highly  permeable  granular  material 
collects  water  from  surrounding,  less  permeable  soil  and 
carries  it  by  gravity  flow  and  slope  to  a  sump  or  outfall, 


where  the  water  can  be  removed  or  wasted.  Today  the  flow 
is  usually  in  slotted  pipes  within  the  granular  trench. 
Since  the  water  can  be  removed  from  the  highly  permeable 
soil  more  quickly  then  it  can  exit  the  less  permeable  soil 
a  difference  in  elevation  of  the  water  levels  of  the  two 
soils  is  created.  This  gradient  serves  to  provide 
continued  flow  as  long  as  the  gradient  exists. 

The  gravity  flow  method  is  obviously  very  easy  to 
implement  and  generally  inexpensive.  It  does  not  always 
work  well  to  dewater  or  drain  a  site  sufficient  to 
accomplish  construction.  It  is  usually  not  effective  to 
attempt  dewatering  by  gravity  flow  when  (Powers,  pg.236): 

a.  Soil  is  a  loose  granular  deposit  without 
plastic  fines. 

b.  The  soil  has  a  high  permeability. 

c.  There  is  a  proximate  source  of  large 
recharge  to  the  water  table  such  as  a  lake  or 
river. 

d.  The  aquifer  is  artesian  or  bounded  under  a 
positive  head  by  an  impermeable  upper  surface. 

e.  The  depth  to  be  dewatered  is  large  so  that 
there  will  be  a  high  gradient  between  the  water 
table  and  the  base  of  the  foundation. 

If  the  soil  to  be  dewatered  does  not  have  the 
limiting  conditions  gravity  dewatering  still  might  not  be 
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advisable.  The  tendency  of  a  soil  to  hold  water  is 
called  storage.  All  soils  will  retain  water  by  capillary 
tension  while  the  apparent  water  table  level  is  drawn 
down.  A  soil  with  a  high  storage  have  a  surprisingly 
large  amount  of  water  held  above  the  water  free 
surface (Powers ,  pg.114).  Removal  of  the  held  water  may 
require  energy  in  the  form  of  pumping. 

For  this  study  the  trench  drain  system  for 
dewatering  beneath  a  foundation  will  be  evaluated.  The 
evaluation  will  be  done  assuming  that  construction  of  the 
foundation  is  complete  and  construction  dewatering  has 
lowered  the  free  water  surface  to  below  the  design  final 
elevation.  The  trench  drain  system  will  maintain  the 
free  water  surface  at  the  final  elevation  and  not  in 
contact  with  the  foundation. 

II  •  Purpose  of  the  Study. 

In  this  study  the  characteristics  of  a  trench  drain 
system  will  be  evaluated  with  the  use  of  the  Aral  Seepage 
Program,  a  finite  element  analysis  program.  Under  a 
range  of  normal  variables  a  trench  drain  system  will  be 
evaluated  to  determine;  (a) the  output  flow  of  the  system 
under  the  different  configurations,  (b)the  range  to  which 
drawdown  of  the  water  table  can  be  expected  and  (c)the 
maximum  rise  of  the  free  water  surface  between  the 
trenches . 
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The  results  of  the  finite  element  analysis  will  be 
grouped  under  dimensional  analysis  to  provide  a  model  for 
prediction  of  the  characteristics.  The  predictive  method 
will  be  effective  within  the  range  of  variables 
considered . 

Finally  a  cost  analysis  of  the  trench  drain  system 
as  compared  to  a  blanket  drain  system.  With  this 
information  a  designer  will  be  able  to  effectively  select 
the  most  cost  efficient  solution  to  the  problem 
considered  if  the  choice  is  between  the  trench  drain  and 
the  blanket  drain  system. 


CHAPTER  2 


PREPARATIONS  FOR  THE  ANALYSIS 


I .  Introduction; 


In  preparation  for  collecting  data  to  analyze  the 
dewatering  capacity  of  a  trench  drain  system  planning  was 
necessary  to  facilitate  the  assimilation  of  results. 
Dimensional  analysis,  used  in  presenting  the  results  of 
this  study  is  briefly  reviewed  in  this  chapter. 

As  with  all  numeric  seepage  calculations  the  area  to 
be  dewatered  is  of  critical  importance.  The  size  of  the 
cone  of  depression  or  the  radius  of  influence  is  the 
single  most  difficult  parameter  to  establish.  Many  rules 
of  thumb,  observations  and  formulas  exist  from  which  the 
radius  of  influence  is  determined  in  practice.  To  predict 
the  rate  of  flow  and  water  table  drawdown  the  radius  of 
influence  must  be  determined.  The  evaluation  of  the 
radius  of  influence  for  use  with  the  Aral  Seepage  Program 
is  reviewed  in  this  chapter. 

The  Aral  Seepage  Program,  a  finite  element  analysis, 
used  in  evaluating  the  trench  drain  system  of  this  study 
was  detailed  by  an  earlier  study  (Pirtle,  Appendix  A)  . 
In  preparing  for  runs  of  this  program  several  items  were 
noted  which  amplify  the  instructions  of  the  user's  manual. 
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II .  Dimensional  Analysis: 

In  this  study  the  modeling  of  a  large  trench  drain 
system  used,  for  dewatering  beneath  a  foundation,  was  done 
using  a  finite  element  analysis  program.  There  are  an 
infinite  number  of  geometries  and  conditions  for  such  a 
system  and  it  is  not  practical  to  run  the  program  for 
each  combination  of  dimensions  and  properties.  The  goal, 
to  provide  a  detailed  summary  of  the  expected  results  for 
any  single  group  of  conditions  from  the  results  of  a  few 
models,  requires  a  range  of  variables  covering  normal 
values.  Results  cannot  be  specific  to  a  few 
configurations  of  the  system. 

Dimensional  analysis  is  a  systematic  grouping  of 
variables  into  a  dimensionless  product.  The  product  can 
represent  a  very  small  model  of  the  system  or  a  full  size 
application.  The  trends  and  results  predicted  from  a 
collation  of  the  dimensionless  numbers  will  be  true  for 
any  combination  of  variables.  From  a  dimensionless  number 
the  quantity  of  flow  or  free  surface  elevation  between 
trenches  can  be  evaluated  for  an  infinite  number  of 
conditions. 

In  arriving  at  the  dimensionless  numbers  for  use  in 
evaluation,  the  variables  must  be  identified  and  cataloged 
by  dimensions  (Langhaar,  Chap  3).  The  variables  are  then 
grouped  into  dimensionless  products.  The  results  of  this 
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study  of  trench  drain  systems  is  presented  in  terms  of 
dimensionless  products.  From  inspection  and 

experimentation  the  dimensionless  numbers  which  provide 
the  most  insight  are  used  in  evaluating  the  models. 

The  variables  in  the  study  of  dewatering  a  foundation 
by  use  of  trench  drains  are  listed: 

-radius  of  influence  (L) ,  units  of  length 
-hydraulic  conductivity,  horizontal  (k^)  ,  units 
of  length/time 

-hydraulic  conductivity,  vertical  (k,,)  ,  units  of 
length/ time 

-trench  spacing,  (S) ,  units  of  length 
-trench  width,  (d) ,  units  of  length 
-free  head  above  trench  bottom,  (h) ,  units  of 
length 

-aquifer  thickness,  (h+H) ,  units  of  length 
-number  of  trench  drains,  (N) ,  no  units 
-slope  of  aquifer,  (p) ,  no  unit 

For  this  study  only  four  trench  drains  were 
considered  so  no  range  of  variables  will  be  available 
from  which  to  reasonably  predict  the  characteristics  of 
systems  with  other  then  four  drains.  The  affect  of  other 
than  four  drains  will  not  be  considered.  With  a 
di\..  i-sionless  product,  later  studies  may  determine  a 
ci  .  elation  between  a  four  trench  system  and  other 
c. onf  •  -^rations. 


The  aquifer  slope  will  often  be  zero  or  near  zero. 
A  multiple  of  zero  will  reduce  any  dimensionless  number 
to  zero  so  the  slope  will  also  not  be  used  in  the 
dimensional  analysis.  Slope  is  a  critical  variable  and 
the  influence  of  aquifer  slope  will  be  shown  by 
individual  results  and  trends.  For  intermediate  results 
interpolation  between  of  the  considered  slopes  will  be 
necessary. 

The  width  of  drainaqe  trenches  considered  for  this 
study  was  restricted  to  one  value.  The  trench  width 
varies  in  practice  and  a  value  of  1.5  feet  provides  the 
minimum  space  necessary  for  placement  of  drainage  pipe 
and  granular  backfill.  If  it  is  more  necessary  to  make 
the  trenches  wider,  the  amount  of  dewatering  will 
increase  (Pirtle,  Chap  4)  and  there  will  be  less 
hydraulic  rise  in  the  water  table  between  trenches.  A 
wider  trench  is  a  more  conservative  approach  and  the 
variable  for  width  of  trench  drains  is  included  as  a 
dimensional  consideration. 

Hydraulic  conductivities,  vertical  and  horizontal, 
are  the  only  variables  with  other  then  a  length  dimension. 
This  complicated  the  forming  of  dimensionless  products. 
The  vertical  and  horizontal  conductivities  had  to  be  in 
the  product,  which  left  five  variables,  all  with  a  single 
length  dimension.  There  are  five  combinations  to  form  a 
dimensionless  product,  each  with  six  alternative 


arrangements.  This  study  was  done  with  the  vertical 
conductivity  (kj  equal  to  the  horizontal  conductivity 
(kj,)  so  the  combinations  for  experimentation  are: 

L,d,S,h  L, d, S , (H+h)  L,d,h, (H+h)  (1) 

L, S , h, (H+h)  S,d,h, (H+h) 

Each  grouping  has  six  different  arrangements, 
yielding  thirty  possible  dimensionless  products.  The 
results  of  the  model  runs  were  prepared  in  the 
dimensionless  form.  The  products  were  compared  and  the 
arrangement  which  provided  an  insight  into  a  relationship 
to  total  flow  and  the  maximum  free  surface  was  selected 
for  use.  In  this  case  the  final  product  is: 

Ld/h (h+H)  (2) 

The  variable  not  included  in  the  product  was  the 
trench  spacing.  In  Chapter  4  trench  spacing  is  considered 
in  evaluating  the  total  flow  and  the  maximum  free  surface 
between  trenches. 

Ill .  Radius  of  Influence: 


The  distance  over  which  the  water  level  is  lowered 
by  dewatering  is  commonly  referred  to  as  the  radius  of 
influence.  The  radius  of  influence  is  the  distance  from 
the  point  or  line  of  water  removal  to  a  point  where  the 
water  table  is  not  affected  by  the  dewatering  operation. 


Obviously  any  model  of  dewatering  or  pumping  operations 
is  dependent  on  the  radius  of  influence. 

The  radius  of  influence  can  be  accurately  determined 
in  the  field  by  a  pumping  test.  A  pumping  test,  using  a 
fully  penetrating  well  in  an  unbounded  aquifer  not  in 
contact  with  a  body  of  water,  will  provide  a  measured 
output  and  drawdown  at  known  radii.  The  drawdown  when 
plotted  against  the  log  of  time  will  allow  for  the 
calculation  of  the  effective  conductivity  of  the  soil. 
Under  the  Dupuit  assumptions  mathematical  models 
((a) Powers,  pg  100)  can  project  a  drawdown  for  various 
distances  and  pump  outputs.  Calculation  of  the  radius  of 
influence  is  also  available  with  this  method.  For  a 
trench  drain  system  the  pumping  test  method  has  several 
problems.  The  trench  drain  is  a  method  of  open  pumping 
or  gravity  flow  and  will  have  a  significantly  lower  output 
then  would  a  pumped  well.  By  its  nature  the  trench  drain 
system  is  longitudinally  aligned  and  would  not  be  modeled 
by  a  well  with  accuracy.  The  trench  drain  does  not  fully 
penetrate  the  aquifer  and  will  offer  a  different  output 
than  a  fully  penetrating  well  (Leonards,  pg  270) . 

It  is  not  unusual  for  the  radius  of  influence  to  vary 
greatly,  often  differing  by  an  order  of  three  to  four 
magnitudes  ((a) Powers,  pg  108).  The  radius  of  influence 
is  affected  by  every  variable  of  the  system.  In  general 
the  radius  will  increase  with  time  and  as  the  drawdown  is 
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increasing.  For  pervious  soils  the  radius  of  influence 
is  greater  then  for  a  soil  of  less  hydraulic  conductivity 
(Leonards,  pg  261) .  In  the  evaluation  of  a  trench  drain 
system,  after  construction  of  the  foundation,  the 
dewatering  has  reached  a  steady  state  condition.  The 
affects  of  time  and  initial  drawdown  will  not  be 
considered  in  this  study. 

Normally  the  radius  of  influence  is  selected  based 
on  knowledge  of  the  area  and  experience  with  design  of 
dewatering  systems.  Several  formulas  exist  to  calculate 
the  radius  of  influence,  usually  considering  the  hydraulic 
conductivity,  thickness  of  aquifer,  free  head  above  point 
of  lowest  drawdown  and  adjustment  factors.  Several 
formulas  are  listed  below: 

Means  of  Calculating  the  Radius  of  Influence 

1.  Q  »  ( .  73+ .  27  (H-h0)  /H)  (kx/2L)  (H2-h2) 

Gravity  flow  to  a  partially  penetrating  slot 
(Leonards,  pg.  271)  with  Q=total  flow, 
H=aquifer  thickness,  h0=elevation  of  water  in 
trench,  x=length  of  trench,  k=  permeability  and 
L=radius  of  influence.  Use  consistent  units. 

2.  L  =  3 . 8h 

Highway  subdrainage  design  (Moulton,  pg  60) 
with  h=drawdown  in  water  level  in  feet, 
L=radius  of  influence  in  feet. 
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3. 


R0=3hk1/2 

Radius  of  influence  as  a  function  of  drawdown 
((a) Powers,  pg.  109  after  Sichart)  with 
R0=radius  of  influence,  h=drawdown.  Use 
consistent  units. 

4.  Q/x  =  K(H2-h2)/L 

Water  table  flow  from  a  line  source  to  a 
drainage  trench  ((a)Powers,  pg.  100).  Q/x  is 
the  total  flow  per  unit  length  of  trench,  k  is 
permeability,  H  =  thickness  of  aquifer, 
h=  trench  elevation,  L  =  radius  of  influence. 

Use  consistent  units. 

5.  G  -  SyLq/KH2 

Predictive  Analysis  of  Groundwater  Inflows  into 
Excavations  (Freeze,  pg.  494). 

G=a  dimensionless  discharge,  Sy=  specific  yield, 
normally  .01  to  .3  (Freeze,  pg.  61), 

K=  hydraulic  conductivity,  H  =  aquifer 
thickness,  q=  rate  of  flow  per  unit  length  of 
trench,  L=  radius  of  influence.  Use  consistent 
units. 

These  formulas  and  others  will  predict  a  radius  of 
influence  for  the  characteristics  of  the  problem.  The 
prediction  of  a  radius  of  influence  for  the  foundation 
trench  drain  system  is  discussed  in  Chapter  4 . 
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The  radius  of  influence  is  an  input  value  in  the 
Aral  Seepage  Program.  The  program  will  calculate  a 
quantity  of  flow  and  a  water  level  free  surface  for  a 
given  radius  of  influence.  The  formulas  from  above  and 
several  others  offer  severe  limitations  to  predicting  the 
radius  of  influence.  The  formulas  are  not  specifically 
for  a  trench  drain  system  and  to  various  degrees  do  not 
account  for  variations  in  horizontal  and  vertical 
permeability ,  partial  penetration  of  the  aquifer  and 
gravity  flow. 

Water  is  produced  from  an  unconfined  aquifer  by 
three  mechanisms:  (1)  expansion  of  the  water  under 
reduced  fluid  pressures,  (2)  compaction  of  the  aquifer 
under  increased  effective  stress  and  (3)  dewatering  of 
the  unconfined  aquifer  (Freeze,  pg  324) .  In  a  gravity 
flow  trench  system  water  is  not  removed  from  the  aquifer 
but,  is  carried  away  after  it  enters  the  trench  (or 
leaves  the  aquifer) .  The  water  will  leave  the  aquifer  at 
a  rate  consistent  with  the  Darcy  principal,  Q  =  kiA. 
Once  steady  state  is  attained  the  area,  A  becomes 
constant  and  so  too  is  the  flow.  The  amount  of  drawdown 
over  the  distance  from  trench  to  original  water  level 
(radius  of  influence)  is  the  hydraulic  gradient,  i.  The 
permeability  is  a  property  of  the  soil  and  is  assumed  to 
remain  constant.  The  gradient  in  a  steady  state  gravity 
flow  system  can  only  change  with  a  rising  (infiltration) 
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or  lowering  (depleted)  water  table.  To  change  the 
gradient  in  a  stationary  water  table  an  outside  action 
must  take  place  to  change  the  fluid  pressures  or  to 
physically  remove  water  from  the  aquifer.  Pumping  would 
be  a  means  of  changing  the  gradient  and  hence  the  radius 
of  influence  and  outflow  of  the  system. 

The  theoretically  correct  radius  of  influence  for 
ideal  conditions  was  determined  by  running  each  trench 
drain  configuration  at  a  series  of  different  radii.  The 
resultant  total  flow  of  the  system  were  plotted  against 
the  trial  radii  of  influence.  Total  flow  is  the 
cumulative  outflow  from  all  four  trenches.  Figures  2-1 
through  2-12  are  the  results  of  the  runs.  The  above 
discussion  shows  the  radius  of  influence  to  be  at  the 
distance  where  no  further  drawdown  occurs.  The  trial 
radii  greater  then  the  true  radius  of  influence  yields  a 
result  at  the  true  free  surface.  Lesser  trial  radii 
results  in  a  forced  solution  with  a  larger  gradient  and 
higher  flow.  For  this  study  the  radius  of  influence  was 
selected  from  the  figures  2-1  through  2-12  as  the 
distance  where  flow  became  constant. 

Establishing  the  distance  where  total  flow  remained 
relatively  constant  required  a  uniform  criterion  to 
define  unchanging  flow.  Arbitrarily,  unchanging  total 
flow  was  defined  as  10%  change  over  a  100  foot  change  in 
radius.  This  definition  was  selected  after  comp^ri '•'■'n  of 


the  results  across  the  thirty  six  different 
configurations  of  the  problem  as  studied.  For  this 
problem,  assuming  a  trench  length  of  200  feet  the  10% 
change  criterion  is  less  then  .05  gpm  total  flow. 

The  radius  of  influence  is  an  extremely  variable 
value,  subject  to  interpretation  and  fine  measurement. 
An  accuracy  of  fifty  feet  was  established  as  a  reasonable 
estimation  for  the  correct  radius  of  influence. 
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Determination  of  radius  of  influence 
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The  Aral  Seepage  Program,  a  finite  element  model 
program  was  used  in  the  evaluation  of  drainage 
characteristics  of  trench  drain  systems.  An  excellent 
program  user's  manual  is  included  in  an  earlier  study 
(Pirtle,  Appendix  A) .  Used  for  confined  and  unconfined 
flow  problems  of  axisymetric  or  two  dimensional  seepage, 
the  program  is  available  at  Georgia  Institute  of 
Technology,  Department  of  Civil  Engineering. 

The  user  inputs  the  problem  geometry,  soil 
characteristics  and  position  of  ground  water.  By 
equilibrium  solutions  to  the  discrete  grid  elements  the 
program  adjusts  the  initially  assumed  free  surface 
through  several  iterations  until  the  level  of  accuracy  is 
satisfied.  When  completed  a  solution  is  presented  to  the 
defined  problem,  including  total  seepage  flow  through 
seepage  faces  and  the  steady  state  free  water  surface. 

The  referenced  user's  manual  includes  detailed 
sample  input  for  an  unconfined  seepage  problem.  The 
evaluated  trench  drain  system  is  of  the  same  nature  and 
so  input  was  extremely  similar.  Appendix  B  of  this  study 
is  a  copy  of  an  input  data  file  for  a  run  of  the  Aral 
Seeoage  Program  for  a  trench  drain.  Figure  3-5  of  Chapter 
3  will  provide  node  and  element  numbers  as  used. 


Comments  on  the  Aral  Seepage  Program  User's  Manual: 

1.  The  user's  manual  does  not  provide  guidance  in 
setting  the  model  sensitivity  on  the  Title  and  Type 
Problem  Card  (card  FERR) .  This  card  sets  the  degree  of 
accuracy  which  will  discontinue  iterations  of  the  free 
surface  calculations.  A  precise  number  is  not  as 
important  as  maintaining  a  minimum  number  of  iterations 
during  the  run.  Recommended  accuracy  will  be  a  value 
such  that  a  minimum  of  three  iterations  are  completed  in 
the  free  surface  calculations.  The  necessary  value  of 
sensitivity  will  become  evident  after  several  runs. 

2.  The  control  card  input  for  NNPC,  total  number  of 
corner  nodal  points,  makes  it  clear  that  intermediate 
nodes  will  be  generated  between  the  listed  corner  nodes. 
What  is  not  clear  is  that  it  is  recommended  that  a  line 
of  elements  several  rows  below  the  free  water  surface, 
should  be  defined  as  corner  nodes.  This  establishment  of 
an  unmoving  row  close  to  the  free  surface  reduces  the 
iteration  time.  Iterations  occur  between  the  intermediate 
nodes  and  the  free  water  surface.  Figure  3-5  is  a  sample 
of  the  finite  element  mesh  used  in  this  study. 

The  Free  Surface  Nodal  Cards  include  the  top  and 
bottom  corner  nodes  of  the  "movable  zone"  described 
above.  Sixteen  corner  nodes,  or  eight  sets  of  free 
surface  and  base,  are  allowed  for  each  card. 
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3.  General  program  notes  not  included  in  the  user's 
manual : 

-The  current  version  of  the  Aral  Seepage 
Program  is  limited  to  thirty  columns  of 
elements . 

-The  total  number  of  seepage  faces  with 
Dirichlet  boundaries  is  not  necessarily  one 
less  then  the  total  number  of  Dirichlet 
Boundary  faces.  In  the  trench  drain  problem 
six  boundary  faces  exist  with  only  four 
drainage  faces.  In  a  problem  done  by  symmetry 
there  would  be  one  less  seepage  face  then  the 
number  of  Dirichlet  boundaries. 

4.  The  output  from  the  Aral  Seepage  Program  is 
extremely  straight  forward  with  only  one  area  of 
confusion.  Appendix  B  of  this  study  includes  a  copy  of 
a  printout  result  from  the  study.  Under  the  "New 
coordinates  of  the  Free  Surface  Line"  there  are  several 
iterations  of  free  surface  corner  nodal  points  and  the  x 
and  y  coordinates.  The  last  iteration  is  the  free 
surface  calculated  within  the  desired  accuracy.  A  plot 
of  the  listed  x  and  y  coordinates  is  a  cross  section  of 
the  free  water  level.  The  maximum  free  surface  elevation 
between  trenches  is  the  maximum  value  between  any  of  the 
four  trench  drains.  In  Chapter  3  the  selection  of  this 
value  is  reviewed.  The  last  page  of  the  output  contains 


the  total  seepage  flow.  The  area  of  confusion  ccmcs  from 
the  presence  of  seepage  output  between  trenches,  across 
element  faces  not  defined  as  seepage  faces.  Those  output 
figures  are  not  correct.  Using  the  nodal  address  numbers 
from  the  bottom  of  the  trench  (Fig.  3-5)  the  actual 
seepage  output  at  each  trench  can  be  obtained.  To  find 
the  total  flow  the  sum  of  trench  flows  is  calculated 
manually  as  the  listed  total  flow  reflects  the  imaginary 
outflow  between  trench  faces. 

V.  Evaluation; 

In  the  analysis  of  the  trench  drain  system 


dimensionless 

products 

will  be 

used 

to 

prepare 

presentations 

of  the 

affect  of 

the 

many 

changing 

variables  upon 

which  the 

solution  is 

dependent . 

The  key 

variable  is  the  radius  of  influence  which  will  be 
determined  by  solving  each  model  for  several  different 
radii  and  selecting  the  correct  value  by  interpolation. 
Finally  in  the  future  use  of  the  Aral  Seepage  Program 
several  notes  can  be  used  to  augment  the  existing  user's 
manual . 
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CHAPTER  3 


INPUT  AND  OUTPUT  OF  THE  TRENCH  DRAIN  EVALUATION 

I .  Introduction: 

The  evaluation  of  different  variables  in  a  foundation 
trench  drain  involves  a  significant  amount  of  input  and 
output  data.  It  is  extremely  important  that  methodical 
presentation  and  collection  of  data  be  practiced.  This 
chapter  details  the  system,  under  which  all  variables  were 
considered.  The  approach  used  insured  that  the  entire 
range  under  consideration  was  evaluated. 

II •  Problem  Geometry: 

In  evaluating  the  dewater*. lg  trench  drain  systems  the 
geometry  must  established.  Among  the  many  variables  in  such 
an  evaluation  are  the  vertical  and  horizontal  hydraulic 
conductivities  of  the  soil ,  number  of  trenches  in  the  system 
and  width  of  the  trench.  To  simplify  the  problem  the  above 
variables  were  held  constant  in  this  study. 

Other  variables  in  an  evaluation  are  spacing  between 
drains,  thickness  of  the  aquifer,  free  head  above  drains  and 
slope  of  the  aquifer.  These  factors  were  varied  in  this  study 
over  ranges  sufficient  to  determine  the  characteristics  of  the 
system  under  as  normally  used. 

Confined  flow  was  not  considered  in  this  study  and  each 
case  was  for  a  water  table  aquifer.  A  confined  aquifer  cannot 
be  dewatered  effectively  with  an  open  pumping  or  sump  method 


((a) Powers,  pg  114).  A  significant  assumption  of  phraetic 
flow,  included  in  the  finite  element  analysis,  provides  the 
radius  of  influence  of  drawdown. 

A  primary  concern  in  dewatering  a  water  table  aquifer  is 
the  storage  depletion  ((c)Powers,  pg.3).  This  study  concerns 
trench  drains  under  steady  state  conditions  following 
construction.  In  steady  state,  the  storage  has  been  depleted 
and  will  not  contribute  additional  flow. 

The  previous  study  (Pirtle)  was  conducted  under  symmetry 
as  a  half  space  problem  with  no  slope.  The  finite  element 
method  in  a  half  space  works  very  well  until  a  sloping  aquifer 
is  considered.  Use  of  symmetry  permits  evaluating  one  half 
of  the  problem,  doubling  the  results  accurately  models  the 
full  problem.  To  do  this  with  a  sloping  aquifer  models  a 
perched  water  table.  Figures  3-1  and  3-2  represents  one  half 
of  a  system  under  symmetry  and  a  full  system  respectively. 

For  this  study  the  affects  of  a  sloping  aquifer  were 
essential,  and  a  full  width  model  was  required.  Figure  3-3 
is  a  representation  of  the  cross  section  of  the  sloping 
aquifer.  In  the  full  width  model  the  free  water  surface  both 
upstream  and  downstream  of  the  foundation  trench  drain  must 
be  evaluated.  Figure  3-4  is  a  cross  section  of  the  problem 
with  the  assigned  variable  names. 

In  the  model  four  trench  drains  were  selected.  Although 
this  models  a  relatively  small  foundation  it  should  yield 
results  which  may  be  used  conservatively  for  larger 


Fig.  3-1  Underdrain  System  in  Half  Space 
(after  Pirtle) 
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Fig.  3-2  Typical  application  of  Trench  Drain 
(after  Pirtle) 
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Fig.  3-3 

Typical  cross  section  of  trench  drain 
system  with  original  aquifer  boundaries 
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Voriables  of  Trench  Orain  Analysis 
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Example  Mesh  Used 


foundations.  Generally  in  a  sloping  aquifer  the  majority  of 
dewatering  and  the  maximum  free  surface  between  trenches  occur 
due  to  the  outboard  drain.  For  a  sloping  aquifer  this  model 
will  closely  model  even  a  large  foundation  as  most  of  the 
removed  water  is  into  the  upstream  trench.  The  number  of 
drains  might  be  a  prime  area  for  future  study. 

The  width  of  the  drainage  trenches  were  held  at  one  and 
one  half  feet  wide.  While  the  earlier  study  (Pirtle,  Chap  4) 
showed  total  flow  to  vary  with  drain  width,  it  is  reasonable 
to  select  a  single,  typical  value.  The  use  of  1.5  feet  is 
based  on  the  approximate  width  of  a  normal  backhoe  bucket. 

The  horizontal  hydraulic  conductivity  (permeability,  kf,) 
was  selected  as  a  constant  0.2  ft2/day/ft  which  is,  7  x  10'5 
cm/sec.  This  is  the  normal  permeability  of  a  silty  sand 
( (b) Cedergren,  pg  34),  common  to  the  Atlanta  area. 

For  the  vertical  hydraulic  conductivity  (Ky)  the  same 
0.2  ft2/day/ft  was  used.  While  a  reasonable  value,  it  does 
not  reflect  the  normal  condition  of  a  greater  horizontal 
hydraulic  conductivity,  which  may  be  greater  than  the  vertical 
by  a  multiple  of  three  to  ten  times  ((a)Powers,  pg.92) .  It 
would  be  an  excellent  study  to  evaluate  the  affects  of 
anisotropy  in  the  trench  drain  system. 

The  primary  area  of  concern  in  this  study  were  the 
affects  of  a  sloping  aquifer.  To  consider  a  complete  range, 
three  slopes  considered.  Initially  0%,  7%  and  15%  slopes  were 
selected.  The  aquifer  thickness  remains  the  same  through  the 
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area  of  consideration  with  a  uniform  aquifer  base  and  water 
table  slope.  Once  the  study  was  begun  it  became  apparent  that 
the  15%  slope  was  too  severe.  The  study  was  modified  with  7% 
the  maximum  slope  considered  and  additional  data  was  collected 
for  a  3.5%  slope. 

Trench  drain  spacing  was  evaluated  using  15  feet,  25  feet 
and  35  feet.  For  the  assumed  four  drain  system  this  allows 
consideration  of  a  foundation  from  60  to  120  feet  wide.  In 
actual  practice  the  trench  drains  are  excavated  between  the 
column  lines  so  normal  spacing  is  that  of  the  column  spacing. 
The  results  for  a  model  with  trench  spacing  of  25  and  35  feet 
are  most  likely  to  be  used. 

An  assumption  for  this  study  was  that  the  aquifer  remain 
of  constant  thickness  across  the  entire  area  of  water  table 
depression.  For  some  models  this  area  is  1000  feet  across. 
It  is  certainly  not  a  normal  soil  condition  for  the  aquifer 
to  be  so  uniform.  In  general  the  area  of  depression  is 
significantly  smaller  then  the  above  value.  Assuming  uniform 
thickness  for  distances  of  300  and  400  feet  is  reasonable. 

To  separate  the  affects  of  free  head  from  thickness  of 
aquifer  the  later  was  considered  as  two  components.  The  first 
component  was  thickness  below  the  drain  or  trench  elevation 
(H)  and  the  second  was  free  head  above  drain  (h) .  The  sum  of 
the  two  components  is  the  total  aquifer  thickness. 

To  bracket  normal  conditions  in  the  evaluation  of  trench 
drains  a  maximum  aquifer  thickness  of  80  feet  was  used.  While 
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arbitrary  this  value  is  reasonable  in  an  aquifer  at  the 
surface.  The  minimum  aquifer  thickness  used  was  30  feet. 
Certainly  smaller  aquifers  exist  at  the  ground  surface  but, 
a  trench  drain  would  probably  not  be  as  effective  as  a  cutoff 
wall  in  protecting  the  foundation  for  such  a  shallow  depth. 

The  last  variable  considered  in  the  study  was  the  free 
head  of  aquifer  above  the  bottom  of  the  trench  drain.  As 
described  above  the  free  head  is  a  component  of  the  aquifer 
thickness.  Free  head  was  held  at  a  maximum  of  20  feet.  As 
the  free  head  increases  the  amount  of  water  removed  must  also 
increase.  Under  an  open  draining  trench  system  the  radius  of 
influence  is  much  larger  and  interference  with  other 
structures  is  a  problem.  As  a  free  head  of  20  feet  is 
approached  alternative  dewatering  methods  may  have  to  be 
considered. 

The  minimum  free  head  considered  was  ten  feet.  This 
value  was  selected  arbitrarily.  The  foundation  trench  system 
would  work  very  well  for  a  smaller  free  head.  In  retrospect 
the  minimum  head  considered  should  have  been  five  feet.  Less 
free  head  and  the  necessary  trench  depth  would  be  so  small 
that  a  blanket  drain  would  be  more  economical.  A  comparison 
of  the  two  drain  systems  will  be  undertaken  on  an  economics 
basis  in  Chapter  4. 
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Ill .  Naming  and  numbering  convention: 


As  a  shorthand  means  of  identifying  the  variables  (Figure 
3-4)  of  this  study  each  is  named  as  follows: 

Radius  of  influence,  L 
Total  flow,  q 
Aquifer  slope,  p 
Trench  spacing,  S 
Trench  elevation,  H 
Free  head  above  trench,  h 
Maximum  free  surface  height,  F 
Hydraulic  conductivity  (vertical) ,  ^ 

Hydraulic  conductivity  (horizontal)  ,  k,, 

As  described  in  Chapter  2  the  determination  of  an 
accurate  radius  of  influence  involved  four  to  six  trial  runs 
of  a  model.  The  result  was  36  geometries  and  over  200 
computer  runs.  In  order  to  identify  the  geometry  of  each 
trial  the  data  files  and  resultants  were  named  from  a  file 
convention.  An  alpha-numeric  name  of  six  characters,  with 


each  one  have  several  values  was  used.  The  files  were  named 
and  the  data  sorted  as  follows: 

A12345.dat,  where  the  .dat  is  a  conventional  suffix 
of  the  operating  system 

First  Character  (A) 

Vertical  hydraulic  conductivity, 

A  =  0.2  ft2/day/ft 
Second  Character  (1) 

Aquifer  slope,  g 
0  =  0% 

3  =  3.5  % 

7  =  7% 

9  =  15  % 

Third  Character  (2^ 

Radius  of  influence,  L 
3  =  100  ft 
5  =  200  ft 

8  =  400  ft 

9  =  500  ft 
0  =  700  ft 

Fourth  Character  (3) 

Trench  spacing,  S 
3  =  15  ft 
5  =  25  ft 
7  =  35  ft 
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Fifth  Character 
Trench  elevation,  H 


2  ■  20  ft 
6  -  60  ft 

Sixth  Character  (5^ 

Free  head,  h 
1  -  10  ft 
2  =  20  ft 

The  above  convention  is  used  in  this  study  in  names  of 
data  files.  If  a  collection  of  data  is  gathered  so  that  each 
group  has  a  constant  trench  elevation,  H  and  free  head,  h  the 
data  groups  will  be  referenced  by  AxxxHh.  In  this  example  if 
the  H  =  20  feet  and  h  =  10  feet  the  data  group  will  be  Axxx21. 
The  data  files  used  in  this  study  are  included  in  Appendix  D. 

IV.  Conduct  of  study: 

The  primary  output  of  the  finite  element  analysis  used 
in  this  study  were  the  total  drainage  outfall  per  lineal  foot 
of  trench  drain  system  and  the  maximum  free  surface  of 
groundwater  between  trench  drains.  As  detailed  in  Chapter  2 
the  study  first  had  to  determine  the  radius  of  influence  for 
each  condition.  From  Figure  3-5  the  elements  of  the  finite 
element  grid  between  can  be  observed.  The  width  between 
trenches  was  divided  into  five  spans.  So  at  four  nodes  the 
iterative  determination  of  the  FEA  program  predicted  an 


ultimate  final  elevation.  This  was  determined  for  every  trial 
radius,  L  and  the  maximum  node  elevation  observed  was 
recorded.  Figures  3-6  through  3-9  are  profiles,  in  the  trench 
area,  of  the  represented  free  surface  for  several  trial  runs. 

It  was  expected  that  the  maximum  free  surface  (F)  would 
occur  near  the  center  of  the  trench  spacing  (S) .  In  order  to 
obtain  an  accurate  maximum  free  surface,  the  nodes  between 
drains  were  not  evenly  spaced.  Instead,  the  distance  between 
nodes  were  greater  near  the  trenches.  For  each  run  the 
maximum  free  surface  selected  was  that  of  the  maximum  nodal 
elevation.  This  is  obviously  not  accurate  in  all  cases.  An 
interpreted  maximum  fre«  surface  (F)  would  increase  many  of 
the  values  used  but,  would  not  reduce  any  value.  The 
anticipated  difference  in  magnitude  of  maximum  free  surface 
(F)  does  not  exceed  0.01  feet  or  0.1  inch.  This  potential 
error  is  not  significant. 

A  review  of  the  maximum  free  surface  value  for  the  many 
different  finite  element  analysis  runs  revealed  that  the 
maximum  free  surface  (F)  always  occurred  between  the  outer 
most  trench  and  the  second  trench.  In  an  aquifer  with  no 
slope  the  conditions  surrounding  the  outer  trench  were  the 
same  for  both  sides.  On  a  sloping  aquifer  the  maximum  free 
surface  (F)  was  at  the  upstream  side  of  the  trench  drain 
system. 

With  the  maximum  free  surface  available  for  every 
geometry  of  problem  and  every  trial  radius  c 2  influence  there 
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FREE  SURFACE  PROFILE 
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remained  the  evaluation  of  the  maximum  free  surface  (F)  at  the 
theoretically  correct  radius  of  influence.  Fortunately  from 
Figures  3-10  through  3-21  it  is  obvious  that  the  maximum  free 
surface  varies  with  trial  radius  of  influence  in  a  consistent 
manner.  As  the  theoretically  correct  radius  of  influence  (L) 
is  known  from  Chapter  2,  the  maximum  free  surface  can  be 
scaled  from  Figures  3-10  through  3-21.  By  convention  the 
interpolation  of  data  to  evaluate  the  radius  of  influence  used 
only  values  which  were  multiples  of  50  feet. 

The  definition  of  a  theoretically  correct  radius  of 
influence  at  the  distance  of  unchanging  total  flow  defines 
the  radius  from  the  known  phraetic  nature  of  the  problem. 
The  unchanging  flow  also  quantifies  the  total  flow  expected 
from  the  system.  For  this  problem  total  flow  was  done  per 
unit  of  width,  making  this  a  two  dimensional  problem. 

The  * otal  flow,  the  maximum  free  surface  and  the  radius 
of  influence  for  each  geometry  of  problem  are  now  available. 
Analysis  of  trench  drain  systems  were  done  with  this  data. 
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Fig,  3-10 

Finding  the  Max.  Free  Surface 
between  trenches 
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Finding  the  Max.  Free  Surface 
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MAXIMUM  FREE  SURFACE 


Fig.  3-13 

Finding  the  Max,  Free  Surface 
between  trenches 
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Fig,  3-14 

Finding  the  Max.  Free  Surface 
between  trenches 
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Fig,  3-15 

Finding  the  Max.  Free  Surface 
between  trenches 
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Finding  the  Max.  Free  Surface 
between  trenches 
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Finding  the  Max.  Free  Surface 
between  trenches 
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Fig.  3-18 

Finding  the  Max.  Free  Surface 
between  trenches 
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Fig.  3-19 

Finding  tho  Max.  Free  Surface 
between  trenches 


MAXIMUM  FREE  SURFACE 

( Slope  7%.  H  =  60',  h  =  10') 


60 


Fig.  3-20 

Finding  the  Max.  Free  Surface 
between  trenches 


MAXIMUM  FREE  SURFACE 


Fig,  3-21 

Finding  the  Max,  Free  Surface 
between  trenches 


V.  Summary ; 

In  summary  there  were  three  slopes  considered,  three 
spacings  between  drains,  two  elevations  of  trench  and  two  free 
heads  above  the  trench.  That  is  thirty  six  geometries  of 
problem,  all  considered  under  one  horizontal  and  vertical 
hydraulic  conductivity.  While  the  physical  dimensions  of  this 
problem  are  adequate  to  cover  the  majority  of  cases  under 
which  trench  drains  would  be  utilized,  the  single  values  of 
k  leave  room  for  further  evaluation.  With  the  use  of 
dimensional  analysis  this  study  might  serve  to  allow 
extrapolation  of  results  for  different  hydraulic 
conductivities  but,  validation  by  further  modeling  is 


recommended . 


CHAPTER  4 


ANALYSIS  OF  TRENCH  DRAIN  DEWATERING 

I .  Introduction; 

Modeling  a  foundation  trench  drain  system  with  the 
a  finite  element  analysis  provides  a  myriad  of  results 
for  many  different  conditions.  For  the  analysis  to  be  of 
use  it  must  be  presented  so  that  future  systems  of  trench 
drains  can  have  capacities  and  limitations  predicted  from 
the  results.  This  chapter  is  the  presentation  of  results 
tor  use  in  prediction  and  selection  of  a  trench  drain 
system. 

II .  Results : 

Chapter  3  details  the  input  and  output  of  this  trench 
drain  analysis.  In  Appendix  A,  Table  A-l  summarizes  the 
different  variables  and  trench  drain  systems  evaluated. 
From  the  data  collected  the  total  flow  of  a  trench  drain 
system  with  four  drains  in  an  isotropic  soil  can  be 
determined  for  a  range  of  design  characteristics  which  are 
reasonable  minimums  and  maximums.  From  this  output 
information  trends  can  be  identified  for  the  response  to 
each  variable. 

The  radius  of  influence  was  found  from  the  measures 
described  in  Chapter  2.  For  each  of  the  thirty  six 
evaluated  geometries  of  trench  drain  system  four  tc  six 


computer  models  were  run,  each  with  a  different  trial 
radius  of  influence.  Figures  2-1  through  2-12  show  the 
results  of  these  trials.  As  expected  each  trial  of  a 
lesser  radius  of  influence  predicted  a  greater  total  flow 
from  the  trench  drain  system.  As  the  trial  radius 
approaches  the  theoretically  correct  radius  of  influence 
the  quantity  of  flow  becomes  unchanging  and  models  the 
gravity  flow  condition  accurately. 

The  spacing  of  trenches  does  not  change  the  total 
flow  of  the  trench  drain  system  in  a  consistent  pattern. 
The  flow  changed  with  other  variables  so  that  the  single 
affects  of  spacing  on  flow  is  indistinguishable. 

It  is  obvious  that  spacing  directly  influences  the 
magnitude  of  the  maximum  free  surface  between  trenches. 
On  Table  A-l  every  combination  of  slope,  trench  elevation 
and  free  head  has  increasing  maximum  free  surface  with 
increasing  spacing.  This  is  consistent  with  expectations. 
Greater  spacing  decreases  the  influence  of  adjacent 
trenches.  Figures  3-6  through  3-9  show  that  the  maximum 
free  surface  in  a  sloping  aquifer  is  strongly  influenced 
by  the  upstream  trench.  The  influence  is  even  more 
strongly  exhibited  with  increasing  space.  A  larger 
spacing  of  trenches  further  isolates  the  upstream  trench 
from  drawdown  influence  from  the  inboard  trenches. 

In  all  cases  the  radius  of  influence  decreased  with 
increased  aquifer  slope.  With  greater  slope,  spacing  of 


the  trenches  exhibited  the  affect  of  decreasing  the  radius 
of  influence.  The  dominant  factor  was  the  slope.  The 
radius  of  influence  was  unchanging  with  spacing  when  the 
slope  was  zero. 

Total  flow  did  increase  with  the  increasing  slope  of 
aquifer  but,  the  magnitude  of  the  change  varied  with  the 
other  variables  so  that  a  single  affect  from  slope  cannot 
be  determined. 

The  maximum  free  surface  increased  consistently  with 
the  aquifer  slope.  The  maximum  free  surface  increased  by 
50%,  for  each  3.5%  incremental  change  of  slope,  in  every 
geometry  considered. 

Impact  of  different  free  heads  above  the  trench  was 
quite  consistent  and  was  normally  coupled  with  the 
elevation  of  the  trench.  For  constant  elevation  of  trench 
the  radius  of  influence  increased  with  free  head.  In  a 
similar  manner  the  total  flow  increased  with  free  head, 
with  the  magnitude  of  change  impacted  by  the  variables 
considered.  The  maximum  free  surface  also  showed  a 
significant  increase  with  free  head  when  the  elevation  of 
the  trench  is  constant. 

For  a  constant  free  head  with  changing  elevation  of 
trench  the  same  trends  in  radius  of  influence,  total  flow 
and  maximum  free  surface  can  be  seen.  The  magnitude  of 
the  change  is  not  as  great  and  the  greater  influence  of 
free  head  can  be  inferred. 
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In  this  study  only  one  drain  width  was  considered 
and  so  the  individual  impact  of  a  changing  width  cannot 
be  determined.  Only  one  vertical  and  horizontal  hydraulic 
conductivity  combination  was  studied  and  all  studies  were 
done  on  a  four  trench  system  so  the  affects  of  these 
single  variables  cannot  be  reported. 

The  many  results  of  this  study  show  the  necessity 
for  a  form  of  dimensional  analysis  in  evaluating  a  model 
affected  by  many  variables.  Some  trends  can  be  predicted 
from  a  couple  of  variables  but,  a  complete  prediction  of 
performance  or  trend  is  not  possible  without  a  dimensional 
analysis. 

Ill .  Results  of  Predictive  Formulas: 

It  is  not  surprising  that  no  single  numeric  method 
exists  for  the  evaluation  of  a  trench  drain  system  of  four 
drains,  in  a  sloping  aquifer.  This  study  attempts  to 
provide  that  solution  but,  do  other,  existing  solutions 
serve  the  need?  In  Chapter  2  several  formulas  were 
referenced  in  the  discussion  of  determining  a  radius  of 
influence.  In  this  section  several  formulas  are 
evaluated.  The  formulas  below  are  presented  with  the 
variables  used  in  this  study  rather  then  as  presented  in 
the  source.  A  comparison  of  the  results  for  the  below 
methods  with  the  results  of  the  finite  element  analysis 
are  in  Table  4-1. 
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a. )  Qp=  (.73  +  .27(h/(h+H))  Kx  (  (h+H)  2-h2) /L 

Qp=total  flow,  h=free  head,  H=trench  elevation,  x=trench 
length,  K=  permeability,  L=  radius  of  influence  (use 
consistent  units) 

The  above  formula  (Leonards,  pg  273)  has  the 
advantages  of  being  for  a  partially  penetrating  slot  in 
an  unconfined  aquifer.  The  disadvantages  are  that  it  has 
no  provisions  for  a  sloped  aquifer,  multiple  trenches, 
trench  width  or  means  for  calculating  the  maximum  free 
surface.  The  radius  of  influence  is  a  value  of  input  to 
this  formula. 

b. )  Figures  36,  61,  68  and  69  (Moulton) 

The  referenced  charts  offer  allowances  for  sloping 
aquifer,  two  drains,  partial  penetration  of  the  aquifer 
and  a  means  of  calculating  a  maximum  free  surface  between 
drains.  Some  disadvantages  of  the  methods  are  the  lack 
of  trench  width,  no  allowance  for  more  then  two  drains  and 
a  half  space  solution  to  the  zero  slope  portion.  The 
radius  of  influence  is  an  input  variable  to  solution.  As 
listed  in  Chapter  3  the  recommended  radius  of  influence 
is  3.8*h.  This  value  yields  a  maximum  radius  of  influence 
of  76  feet  for  this  study.  This  value  was  much  less  then 
used  throughout  the  study  and  yielded  results  far  greater 
then  predicted  from  the  finite  element  analysis.  In  Table 
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4-1  the  radius  of  influence  used  was  that  from  the  finite 


element  analysis  which  yielded  results  much  closer  to 
those  of  the  other  calculations. 

c.)  Q/x  =  K((H+h)2  -  h2)/L  ((a)  Powers,  pg  100) 

Q=  total  flow,  x=  trench  length,  K=  permeability, 
H=  trench  elevation,  h=  free  head,  L=  radius  of  influence 
(  Use  consistent  units) 

This  method  is  a  traditional  trench  solution.  It 
does  not  make  allowance  for  partial  penetration  of  the 
aquifer,  multiple  trenches,  trench  width  or  aquifer  slope. 


d.) 


G  =  S  Lq/KH2  and  figures  18(b)  and  (c) 


(Freeze,  pg  261)  G=  dimensionless  product,  Sy=  Specific 
Yield,  L=  radius  of  influence,  q=  total  flow  per  unit 
width,  K=  permeability,  H=  aquifer  thickness  (use 
consistent  units) 

This  solution  is  for  a  partially  penetrating  trench. 
Once  again  this  method  makes  no  provision  for  multiple 
trenches,  trench  width,  sloping  aquifer.  In  addition  the 
Specific  Yield  is  used  adding  another  variable,  this  one 
normally  varying  between  .01  and  .3  (Freeze,  pg  61). 

.  Each  of  the  above  methods  have  several  disadvantages 
and  were  not  developed  for  use  in  the  evaluation  of  a 
trench  drain  system.  Assuming  that  the  finite  element 
analysis  is  the  most  accurate  of  the  means  of  evaluation 
tried  it  is  obvious  in  Table  4-1  that  no  single  formula 
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TABLE  4-1  Summary  of  Predictive  Methods 


Summary 

A 

B 

<2 

n 

SI  ope 

0 

3.5 

3.5 

7 

Spacing  (ft) 

35 

15 

35 

25 

Trench  El  (ft) 

60 

20 

60 

20 

Free  Head  <  f  t ) 

10 

20 

20 

10 

(Mansur ,  pg  273) 

total  f  1  ow  (■ft'' 2/da y) 

0.5 

0.46 

0.81 

0.55 

(Moulton,  pg  120  &  130) 

total  flow  (ftA2/day) 

0.49 

0.22 

0.73 

0.42 

(Powers,  pg  100) 

total  flow  (ftA2/day) 

0 . 66 

0 . 53 

1 .02 

0 . 66 

(Freeze,  pg  261) 
total  flow  (ftA2/day) 

5.88 

0.64 

5.12 

O 

Am 

Aral  Seepage  Program 

total  flow  (ftA2/day) 

0.14 

0.19 

0 . 25 

0 . 27 

Radius  of  Influ  (ft) 

500 

500 

500 

150 

Trial  Radius  200'  ^ 

Total  Flow  (ft  /day) 

•  3 

.45 

.54 

.26 

approximates  the  solution  to  the  problem  very  closely. 
Of  the  evaluated  methods  the  Moulton  method,  using  the 
radius  of  influence  as  determined  by  the  finite  element 
analysis  was  the  most  accurate.  It  is  worth  noting  that 
with  no  other  means  of  selecting  a  radius  of  influence  the 
value  determined  in  this  study  was  used  for  each  method. 
This  approximation  would  certainly  affect  the  results. 
Appendix  E  is  a  discussion  of  the  affect  of  arbitrary 
selection  of  radius  on  total  flow  and  free  surface  height. 
IV.  Prediction  of  Trench  Drain  Performance: 

The  prediction  of  total  flow  and  maximum  free  surface 
is  possible  through  the  use  of  results  presented  in  a 
dimensionless  form.  As  a  trial  some  runs  were  done  to 
explore  the  affects  of  anisotropy  and  number  of  trenches. 
These  results  are  discussed  below. 

Any  prediction  of  flow  characteristics  requires  an 
evaluation  of  the  radius  of  influence.  Figures  4-1 
through  4-3  are  useful  for  predicting  the  radius  of 
influence  from  the  trench  elevation,  free  head,  slope  of 
aquifer  and  the  trench  spacing.  These  figures  are  only 
appropriate  for  an  isotropic  condition  in  a  four  trench 
system  with  a  trench  width  of  1.5'.  Use  of  the  figures 
will  frequently  require  several  interpolations. 
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Finding  the  radius  of  influence 


To  use  figures  4-1  through  4-3  the  following  example 
will  be  useful. 

Example  4-1:  Find  the  radius  of  influence  for  a  four 
trench  system  with  trench  width  of  1.5'  and  = 
1  with  K  =  .2  ft/day.  The  trench  elevation  is  40 
feet,  the  free  head  is  15  feet,  the  trench  spacing 
is  30  feet  and  the  aquifer  slope  is  4%. 

a.  Figure  4-4  is  annotated  from 
Figure  4-2.  Follow  this  figure 
through  the  following  steps. 

b.  At  slope  =  4%  on  Figure  4-2 
interpolate  between  the  lines  labeled 
2  0/10  and  20/20  as  well  as  between 
60/10  and  60/20.  These  are  the 
points  for  trench  elevations  of  20 
and  60  feet  with  a  free  head  of  15 
feet.  The  radii  of  influence  for 
these  two  points  are  330  feet  and  456 
feet. 

c.  Again  at  slope  =  4%  interpolate 
between  the  above  two  points  to  find 
the  radius  of  influence  for  a  trench 
elevation  of  40  feet  and  a  free  head 
of  15  feet  with  a  spacing  of  25  feet. 

This  radius  of  influence  is  394  feet. 
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Radius  of  Influence 

Four  trench  drains.  Slope  0  to  7% 


nding  the  radus  of  influence 
>r  example  4  -1 


d.  Follow  the  above  procedure  again 
on  Figure  4-3  to  get  a  radius  of 
influence  of  344  feet  when  spacing  is 
35  feet. 

e.  For  a  spacing  of  30  feet 
interpolate  between  the  results  of 
steps  c.  and  d.  above  to  get  a  final 
predicted  radius  of  369  feet. 

As  an  alternate  method  a  rough  selection 
of  350  or  400  feet  will  not  make  a  significant 
difference  in  the  following  calculation  method. 

As  discussed  in  Chapter  2  the  radius  of 
influence  is  an  extremely  variable  value,  often 
approximated  only  within  an  order  of  magnitude. 

One  suggested  method  (Moulton,  pg  66)  has  the 
radius  of  influence  estimated  as  3.8  times  the 
free  head.  In  this  example  that  would  predict 
a  radius  of  only  53  feet,  one  seventh  of  that 
from  this  study  but,  within  one  order  of 
magnitude. 

In  Chapter  2  the  derivation  of  the  dimensionless 
product  used  in  this  study  was  reviewed.  The  product 
which  best  presented  the  behavior  of  total  flow  and 
maximum  free  surface  with  the  changing  variables  of  the 
problem  geometry  was: 
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Ld/h (h+H)  (2) 

Where  L  =  radius  of  influence,  d  =  trench  width,  h  =  free 
head  and  H  =  trench  elevation  as  shown  on  Figure  3-4. 
Using  consistent  units  for  the  variables  yields  a 
dimensionless  product. 

Figures  4-5  and  4-6  are  total  flow  and  maximum  free 
surface,  respectively  against  the  dimensionless  product 
for  the  considered  geometries.  In  both  plots  an 
acceptable  curve  fit  is  possible.  The  points  are  aligned 
with  the  largest  aquifer  slope  to  the  left  and  the  zero 
sloped  aquifer  to  the  right.  This  observation  is  for 
interest  only  as  the  dimensionless  product  provides 
sufficient  accuracy  with  the  sloping  aquifer  influence 
represented  by  the  values  of  total  flow,  maximum  free 
surface  a;  ''  radius  of  influence.  From  the  previous 
example : 

a.  Calculate  Ld/h (h+H).  For  this  problem 

using  the  above  L  =  369  feet, 

Ld/h ( h+H)  =  369 ' *1. 5 '/I5  '  (15 ’+40  '  )  =  .671 

b.  From  Figure  4-5  total  flow  q  =  .25  ft?/day 

c.  From  Figure  4-6  maximum  free  surface 

F  =  .9  ft 

As  a  check  of  the  accuracy  of  this  method  to  predict 
the  output  of  the  finite  element  analysis  six  sample  runs 
were  conducted.  The  values  of  the  variables  and  the 
results  of  the  finite  element  analysis  are  shown  in  Table 
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Predicting  flow  from  a 
dirfiensionless  product 
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Rinding  the  maximum  free  surface 
freon  a  dimensionless  product 


A-2.  For  this  analysis  three  trial  radii  of  influence 
were  used  for  each  sample  run.  The  theoretically  correct 
radius  of  influence,  total  flow  and  maximum  free  surface 
were  taken  from  plots  similar  to  Figures  2-1  and  3-10. 

The  determination  of  the  theoretically  correct 
radius  of  influence  is  done  under  an  arbitrary  standard 
to  an  accuracy  of  fifty  feet.  From  Table  A-2  the 
prediction  of  total  flow  and  maximum  free  surface  agree 
with  the  results  of  the  finite  element  sample  runs.  This 
is  excellent  accuracy  in  selecting  the  radius  of 
influence.  This  method  is  acceptable  for  an  isotropic 
soil  with  a  four  trench  system. 

Figures  4-5  and  4-6  were  prepared  assuming  an 
isotropic  condition.  As  a  matter  of  interest  several 
finite  element  runs  were  done  modeling  a  four  drain 
system  with  various  conditions  of  anisotropy.  Normally 
soils  are  anisotropic  with  KyK,,  between  .1  and  .33.  The 
results  of  these  sample  runs  are  shown  on  Figures  4-7  and 
4-8,  where  the  affects  of  differing  hyuraulic 
conductivity  and  anisotrophic  soil  are  explored. 

From  Figure  4-7  it  appears  that  different  lines, 
nearly  parallel ,  exist  for  a  four  dra-  .stems  with  the 
line  of  largest  total  flow  having  the  highest  hydraulic 
conductivity  and  the  line  of  lowest  flow  with  the  lowest 
hydraulic  conductivity.  The  lines  of  differing 
conductivity  may  converge  to  the  left  side  of  the  plot. 
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MAXIMUM  FREE  SURFACE 


ipy  and  maximum  free  surface 
dimensionless  product 


The  degree  of  anisotropy  appears  to  have  no  consistent 
pattern. 


Figure  4-8  shows  the  opposite  influence.  The  degree 
of  anisotropy,  represented  by  Kyi^,  would  have  parallel 
lines  with  decreasing  maximum  free  surface  as  the  degree 
of  anisotropy  decreases.  The  lines  of  this  plot  may  also 
converge  to  the  left  as  Ld/h(h+H)  decreases. 

Another  area  of  interest  is  the  affect  of  a 
differing  number  of  trench  drains.  As  a  means  of 
predicting  the  ultimate  affect  of  systems  with  more  or 
less  drains  then  the  four  evaluated.  Figures  4-9  and 
4-10  show  the  results  of  three  runs  of  a  three  trench 


system.  While  not  conclusive  it  is  likely  that  the  line 


for  a  three  trench  system  is  just  below  the  four  trench 


line.  With  further  study  it  might  be  possible  to  extend 
the  results  of  a  four  drain  analysis  based  on  the 


paralleling  trends  of  lines  for  the  different  trench 


configurations . 
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MAXIMUM  FREE  SURFACE 
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3  &  4  Trenches  for  maximum  free  surface 
frotn  a  dimensionless  product 


This  study  has  developed  an  effective  method  for 
determining  the  total  flow  and  maximum  free  height  for  a 
trench  drain  system,  knowing  the  geometry  of  the  problem. 
It  appears  probable  that  anisotropic  conditions, 
different  hydraulic  conductivities  and  trench  systems 
with  other  then  four  drains  can  be  included  into  this 
type  of  evaluation  with  some  further  research. 

V.  Comparison  of  Trench  Drain  with  Blanket  Drain: 

An  alternative  to  the  trench  drain  system  for 
dewatering  beneath  a  foundation  is  the  blanket  drain.  A 
blanket  drain  system  is  similar  to  the  trench  drain 
system  as  it  is  a  mass  of  open  graded  stone  or  gravel 
isolated  from  the  soil  by  a  geosynthetic  filter  material 
and  has  slotted  collector  pipes  to  carry  water  to  a  sump 
for  removal.  It  is  different  in  geometry.  The  blanket 
drain  is  a  single  rectangular  cross  section  which 
dewaters  the  surrounding  soil  to  the  depth  of  the  bottom 
of  the  section.  The  trench  drain  uses  less  gravel  by 
extending  the  dewatering  depth  down  with  narrow  trenches 
extending  from  a  thin  base  of  gravel  across  the  bottom  of 
the  foundation  (Figure  3-4) . 

The  blanket  drain  does  an  excellent  job  of 
dewatering  beneath  a  foundation  and  has  the  advantage  of 
easier  construction.  The  biggest  difference  between  the 
two  systems  in  a  practical  application  is  the  cost.  The 
trench  system  requires  more  detailed  labor  and  excavation 
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while  the  blanket  drain  requires  more  gravel.  To  do 
a  cost  comparison  of  the  two  methods  a  standard 
foundation  size  was  necessary  to  determine  the 
excavation,  placement  and  compaction  productivity  rates. 
In  this  study  a  foundation  of  100  feet  by  200  feet  was 
used,  only  for  the  establishment  of  unit  costs  of  labor. 

The  evaluation  is  aone  based  on  a  set  slab  elevation 
and  a  predetermined  base  of  trench  or  blanket  drain. 
That  is  the  blanket  drain  will  have  to  be  thick  enough  to 
dewater  to  the  depth  of  the  trench  drain  without  lowering 
the  foundation.  A  basic  cost  to  this  evaluation  was  the 
cost  of  gravel  (#57  stone)  which  goes  for  $8. 30/ton  in 
the  Atlanta  area.  Labor  costs  used  were  $20/hr  for  an 
equipment  operator  and  $16/hr  for  a  laborer  (both  rates 
including  overhead  and  fringes) .  Productivity  values 
came  from  the  Means  estimating  guide  for  1989.  The  cost 
of  equipment  rentals  were  approximated  from  the  same 
guide.  The  cost  of  the  geosynthetic  filter  material  was 
estimated  as  $0. 70/yd2  from  discussions  with  suppliers  and 
engineers  in  the  Atlanta  area. 

The  cost  data  was  input  into  a  spreadsheet  with  a 
variety  of  possible  trench  drain  configurations.  The 
costs  were  calculate^  for  both  the  trench  drain  and  an 
equivalent  blanket  drain  for  several  different  systems 
each  with  more  drains.  Table  4-2  is  a  copy  of  one  such 
trial.  With  the  cost  of  gravel  so  much  higher  then  the 
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filter  material  the  ratio  of  trench  drain  cost  to  blanket 
drain  cost  was  less  then  one  in  all  cases. 

With  a  blanket  drain  at  a  shallow  elevation  the  cost 
comparison  is  more  likely  to  recommend  a  blanket  drain 
system.  Figure  4-11  is  a  summary  of  such  a  comparison. 
To  use  a  thickened  blanket  drain  will  sacrifice  drawdown 
capacity  in  exchange  for  some  construction  ease  and  cost 
benefits.  The  decision  will  depend  on  the  application. 
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TABLE  4-2 


Trench  and  Blanket  Drain  Costs 


Ex  c  a^  a  t  i  on 
Trench 

B 1 anko  t 


unit  costs 

0.4  $/cu .i t 

0.08  $/cu  .  -f  t 


Spacing  Trench  Trench 
Width  Depth 


Material  costs  (  i  i 
Geosyn  0.08 

Gravel  Tr  0.45 
Gravel  61  0.9 

per  unit 
Number  o-fTrench 
Trenches  Total  $ 


1  5 

1  .5 

1 

30 

38 . 25 

1  5 

1  .5 

4 

30 

153 

15 

2.5 

1 

30 

63.75 

15 

2.5 

4 

30 

255 

25 

1  .5 

1 

30 

38.25 

25 

1  .5 

4 

30 

1  53 

■-a 

2.5 

1 

30 

63 . 75 

25 

2.5 

4 

30 

255 

35 

1  .5 

1 

■Sf  0 

38  •  25 

35 

1 .5 

4 

30 

153 

35 

2.5 

1 

oo 

w 

63.75 

35 

2.5 

4 

30 

255 

place) 
$/sq  .  1 1 

$/cu  .  1 1 
$/cu  .  +■  t 

1  n  g  t  h 
B1 anke  t 
Total  % 

508 . 96 
1 920 . 84 
540 . 76 

2040.64 
818. 38 

3080 . 64 
848 . 1 6 

3200 . 64 
1123.76 

4240.64 
1155.56 

4360 . 64 
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Spacing  for  ecconomic  use  of  Trench 
Drdin,  dewatering  to  greater  depth 


VI .  Conclusions : 


In  this  chapter  a  method  for  predicting  the  results 
of  a  finite  element  analysis  of  a  trench  drain  system  for 
isotropic  soil  and  four  drains  is  presented.  A  rough 
approximation  has  been  made  to  show  the  possible  impact 
of  anisotropy  and  different  trench  configurations. 
Empirical  data  is  not  available  to  validate  or  modify  the 
predictions  of  the  Aral  Seepage  Program  for  a  trench 
drain  system.  Compared  with  other  slot  drain  solutions 
it  appears  that  the  finite  element  solution  predicts 
lesser  flows  and  is  so  less  conservative. 

The  use  of  a  trench  drain  system  can  be  justified 
economically  over  a  blanket  drain  system  in  some 
configurations.  The  brief  evaluation  of  this  chapter 
provides  a  simple  comparison  method  for  the  relative  cost 
differences  of  the  two  systems. 
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CHAPTER  5 


CONCLUSIONS 


I .  Introduction: 

In  this  study  the  trench  drain  system  was  evaluated  for  thirty 
six  configurations  using  a  finite  element  analysis  program.  The 
results  of  the  analysis  were  used  under  dimensional  analysis  to 
form  a  method  of  predicting  the  results  of  a  similar  analysis  for 
any  configuration  of  a  trench  drain  system  within  the  limits  of  the 
study.  The  final  portion  of  the  study  was  to  evaluated  the  cost 
of  a  trench  drain  system  as  compared  with  the  blanket  drain  system 
for  use  in  a  drain  beneath  a  foundation. 

II .  Trench  Drain  System  Analysis: 

The  results  of  the  finite  element  analysis  could  not  be 
validated  by  test  or  field  data  nor,  did  the  formulas  available 
show  close  agreement  with  the  results  of  the  analysis.  It  would 
be  useful  if  the  actual  outflow  of  a  trench  drain  system  were 
measured  and  the  radius  of  influence  of  drawdown  determined. 

This  study  was  limited  to  isotropic  conditions  and  a  four 
drain  system.  The  affect  of  other  conditions  was  briefly  explored 
but,  a  more  complete  study  is  necessary  before  the  results  are 
extrapolated  to  conditions  of  many  drains  or  variable  soils. 
Initial  expectations  are  that  results  of  Figures  4-7  through  4-10 
will  be  consistent  with  results  of  a  more  complete  run  of 


variables.  Whether  more  then  four  drains  have  significantly 
different  shape  then  the  four  and  three  drain  curves  will  be  most 
important  when  the  cost  comparison  of  trench  and  blanket  drains  is 
considered  from  Figure  4-11. 

The  one  recurring  constant  in  the  reference  material  on 
dewatering  is  the  variable  nature  of  the  radius  of  influence.  For 
the  narrow  conditions  evaluated  in  this  study  reasonable 
estimations  of  the  radius  of  influence  are  possible  from  Figures 
4-1  through  4-3.  The  earlier  recommendations  of  field  observations 
and  computer  modeling  with  other  ranges  of  variables  would  expand 
the  information  available  for  predicting  the  radius  of  influence. 

The  permeability  used  in  this  study,  0.2  ft/day  is  consistent 
with  the  soil  in  the  Atlanta  area.  If  the  trench  drain  system  is 
to  be  used  in  a  greatly  more  permeable  soil  the  results  may  be 
significantly  different.  It  would  not  be  recommended  to  use  a 
trench  drain  in  a  very  porous  sand  or  gravel  just  as  an  open  sump 
is  not  satisfactory  to  dewater  a  site  of  that  composition. 
Establishing  the  upper  and  lower  bounds  of  permeability  for 
practical  use  of  trench  drains  is  a  most  useful  recommendation  for 
further  study. 

Ill .  Predictive  Method: 

The  predictive  method  using  Figures  4-1  through  4-6  work  well 
in  predicting  the  results  of  the  finite  element  analysis.  The 
forms  are  easy  to  use  and  adapt  readily  to  drain  configurations 
within  the  range  evaluated.  The  range  of  variables  were  selected 


to  match  the  conditions  under  which  a  trench  drain  system  is 
normally  used.  Extrapolation  of  the  results  beyond  these  values 
may  be  of  limited  practicality  and  questionable  accuracy. 


The  prediction  of  the  maximum  free  surface  at  a  low 
dimensionless  product  (Ld/h(h+H))  should  be  treated  with  great 
care.  When  the  product  is  less  then  0.5  the  curve  rises  rapidly 
and  the  maximum  free  surface  exceeds  1.25  feet  only  in  this  range. 
The  maximum  free  surface  is  of  concern  only  when  the  water  level 
approaches  the  shallow  blanket  immediately  beneath  the  slab, 
approximately  one  foot  for  most  trench  systems.  The  curve  for  the 
maximum  free  surface  was  roughly  aligned  with  the  7%  slope  values 
at  the  left  end  and  the  0%  slope  values  on  the  right.  For  large 
aquifer  slopes  the  curve  may  approach  vertical  and  might  be 
sufficient  reason  to  avoid  the  use  of  a  trench  drain  system  in  such 
an  aquifer. 

IV.  Cost  Comparison: 

If  the  necessary  depth  of  dewatering  is  to  the  depth  of  the 
trench  drain  system  a  blanket  drain  cannot  compete  on  an  economic 
basis.  If  a  very  shallow  level  of  dewatering  is  required  beneath 
the  foundation  it  would  not  be  practical  to  install  ~  t -ench  drain 
system.  Trenches  of  under  one  foot  depth  would  requ.<  individual 
consideration.  By  the  results  of  Figure  4-11  it  is  apparent  that 
even  sacrificing  the  depth  of  dewatering  is  not  sufficient  to  turn 
the  advantage  to  blanket  drains  in  all  cases. 


* 


a  •* 


V.  Evaluation: 


The  results  of  this  study  have  satisfied  the  stated 
intentions.  Results  and  input  were  reviewed  with  care  and  show 
consistent  and  reasonable  trends.  Appendix  C  is  a  brief 
description  of  the  method  of  compilation  of  data,  preparation  of 
data  files  and  execution  of  program  commands.  Also  included  in 
Appendix  D  is  a  storage  disk  for  an  IBM  compatible  computer  with 
the  spreadsheets  of  results  and  the  data  files  used. 
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APPENDIX  A 

SUMMARY  OF  PROGRAM  COMPUTER  RUNS 


APPENDIX  A 


FINITE  ELEMENT  COMPUTER  RUN  SUMMARY 


8pac i og 

Tr*neh 

El *vat i on 

Fr** 

H»ad 

Radius  oi 
In* 1 u*nc* 

Total 

Flow 

Max  Fr** 
Surlac* 

Numb*r  o* 
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81  op*  ■  0 V. 

15 

20 

10 

400 

0.11 

0.18 
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15 

20 

20 

550 

0.2 

0.33 
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15 
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500 

0.17 

0.35 
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15 
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0.29 
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400 
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0.35 
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60 
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0.15 
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5 

35 

20 

10 
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5 

35 

20 

20 
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0.38 

5 

35 

60 

10 

500 

0.14 
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5 
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60 

20 
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5 

Slop*  ■  3.5 X 
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20 
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0.14 

0.38 

4 

15 

20 

20 

500 

0.19 

0.52 

4 

15 

60 

10 

450 

0.15 

0.53 

4 

15 

60 

20 

550 

0.26 

0.93 

4 

25 

20 

10 

250 

0.14 

0.5 

4 

25 

20 

20 

450 

0.18 

0.69 

4 

25 

60 

10 

400 

0.16 

0.79 

4 

25 

60 

20 

550 

0.24 

1.27 

4 

35 

20 

10 

200 

0.15 

0.66 

4 

35 

20 

20 

400 

0.18 

0.82 

4 

35 

60 

10 

350 

0.18 

1  .02 

4 

35 

60 

20 

500 

0.25 

1.35 

4 

SI  op*  •  T/. 

IS 

20 

10 

200 

0.23 

0.6 

5 

15 

20 

20 

350 

0.3 

0.82 

5 

15 

60 

10 

300 

0.23 

0.819 

5 

IS 

60 

20 

500 

0.29 

1  .016 

5 

25 

20 

10 

150 

0.27 

0.85 

5 

25 

20 

20 

250 

0.32 

1.17 

5 

25 

60 

10 

200 

0.28 

1.27 

5 

25 

60 

20 

450 

0.31 

1  .4 

5 

35 

20 

10 

100 

0.28 

1.14 

5 

35 

20 

20 

200 

0.33 

1.432 

5 

35 

60 

10 

150 

0.33 

1  .69 

5 

35 

60 

20 

400 

0.34 

1  .85 

5 

Table  A-l  Summary  of  Finite  Element  Runs 
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SAMPLE  INPUT  AND  OUTPUT 
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SAMPLE  INPUT  DATA  FILE 


“A35521 “ 

2,0.00001 

90,29,0,30,20,6,4,0 

I , -190,6.65, 1 
7,-190,22.65,1 

I I , -190 ,36.65,0 
18,-114,3.99,1 
24,-1 14, 19.99, 1 
28,-114,32.99,0 
35,-43.7,1 .5295,1 
41 ,-43.7,17.5295,1 
45,-43.7,29.5295,0 
52,-19,0 .665, 1 
58,-19,16.665,1 
62,-19,27.665,0 
69,0,0,1 
75,0,16,1 
79,0,25,0 
86,10,-0.35,1 

92.10.16.1 
96,10,20,0 

103.11 .5, -0.4025,1 

109.11.5.16.1 
113,11.5,20,0 

120.16.5, -0.5775,1 

126. 16.5. 16. 1 
130,16.5,22,0 
137,23,-0.805,1 

143.23.1 6.1 
147,23,22,0 
154,25,-0.875,1 

160 .25. 16. 1 
164,25,22,0 

171  ,31 .5,-1 .1025,1 

177.31 .5.16.1 
131,31 .5,22,0 

183.36.5, -1 .2775,1 

194.36.5. 16. 1 
198,36.5,20,0 
205,38,-1 .33,1 

211.38.16.1 
215,38,20,0 
222,43,-1 .505,1 

228.43.16.1 
232,43,22,0 

239.49.5, -1 .7325,1 

245.49.5.16.1 
249,49.5,22,0 

256.51 .5, -1 .8025,1 

262.51 .5.16.1 
266,51 .5,22,0 
273,58,-2.03,1 

279.58.16.1 
283,58,22,0 
290,63,-2.205,1 
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296.63. 16. 1 
300,63,20,0 

307.64.5, -2.25/5, 1 

313.64.5. 16. 1 
317,64.5,20,0 

324.69.5, -2.4325,1 

330 .69.5. 16. 1 
334,69.5,22,0 
341 ,76,-2.66, 1 

347.76. 16. 1 
351 ,76,22,0 
358,78,-2.73,1 

364.78.16.1 
368,78,22,0 

375.84.5, -2.9575, 1 

381 .84.5.16.1 
385,84.5,22,0 

392.89.5, -3.1325,1 

398.89.5.16.1 
402,89.5,20,0 
409,91 ,-3.185,1 

415.91 .16.1 
419,91 ,20,0 
426,101 ,-3.535,1 

432.101 .15.65.1 
436,101 ,21.465,0 
443,120,-4.2,1 

449.120.14.985.1 
453,120,21 .8,0 

460 , 144.7,-5.0645, 1 

466. 144.7. 14. 1205. 1 
470 , 144.7,22.9355,0 
477,215,-7.525,1 

483.215.1 1 .66. 1 
487,215,21 .475,0 
494,291 ,-10.185,1 

500.291 .9.1 
504,291 ,19.815,0 

1,12,18,19,20,13,3,2,4,0.2,0.2,0 
1 8 , 29 , 35 , 36 , 37 , 30 , 20 , 1 9 , 4 , 0 . 2 , 0 . 2 , 0 
35,46,52,53,54,47,37,36,4,0 .2,0 .2,0 
52,63,69,70,71 ,64,54,53,4,0.2,0.2,0 
69,80 ,86,87,88,81 ,71 ,70 ,4,0 .2,0 .2,0 
86,97, 103, 104, 105,98,88,87,4,0 .2,0 .2,0 
103,114,120,121,122,115,105,104,4,0.2,0.2,0 
120,131 ,137,138,1 3? ,132,1 22 ,121 ,4, 0.2,0. 2,0 
1  37 , 1 48 , 1 54 , 1 55 , 1 56 , 1 49 , 1 39 , 1 38 , 4 , 0 . 2 , 0 . 2 , 0 
154, 165, 171 , 172, 173,166, 156, 155,4,0.2,0.2,0 
171,182,188,189,190,183,173,172,4,0.2,0.2,0 
188,199,205,206,207,200 ,190,189,4,0.2,0.2,0 
205,216, 222 , 223 , 224 ,217,207,206,4,0.2,0.2,0 
222 , 233 , 239 ,240,241, 234 , 224 , 223 ,4, 0.2, 0.2,0 
239,250,256,257,258,251 ,241,240,4,0.2,0.2,0 
256 , 267 , 273 , 274 , 275 , 268 , 258 , 257 ,4, 0.2, 0.2,0 
273,284,290,291 , 292 , 285 , 275 , 274 , 4 , 0 . 2 , 0 . 2 , 0 
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I 


# 


I 


290,301 ,307,308,309,302,292,291 ,4, 0.2, 0.2,0 
307,318, 324 , 325 , 326 ,319,309,308,4,0.2,0.2,0 
324 , 335 , 341 , 342 , 343 , 336 , 326 , 325 ,4, 0.2, 0.2,0 
34 1 , 352 , 358 , 359 , 360 , 353 , 343 , 342 ,4, 0.2, 0.2,0 
358 , 369 , 375 , 376 , 377 , 370 , 360 , 359 ,4, 0.2, 0.2,0 
375 , 386 , 392 , 393 , 394 , 387 ,377, 376 ,4,0. 2, 0.2,0 
392,403,409,410,41 1 ,404,394,393,4,0.2,0.2,0 
409,420,426,427,428,421 ,411,410,4,0.2,0.2,0 
426 , 437 , 443 , 444 , 445 , 438 , 428 , 427 ,4, 0.2, 0.2,0 
443 , 454 , 460 , 46 1 , 462 , 455 , 445 , 444 ,4, 0.2, 0.2,0 
460,471 ,477,478,479,472,462,461 ,4, 0.2, 0.2,0 
477 , 488 , 494 , 495 , 496 , 489 , 479 , 478 ,4, 0.2, 0.2,0 

1 .36.65.5.2. 1 

96.20 . 1 . 17.6 

198.20.1 .17.6 

300.20 .1.17.6 

402.20.1 .17.6 

494.19.815.5.2.1 

1 1 , 7 , 28 , 24 , 45 , 41 , 62 , 58 , 79 , 75 , 96 , 92 , 1 1 3 , 1 09 , 1 30 , 1 26 
1 47 , 1 43 , 1 64 , 1 60 , 1 81 , 1 77 , 1 98 , 1 94 , 21 5 , 21 1 , 232 , 228 , 249 , 245 , 266 , 262 
283 , 279 , 300 , 296 ,317,313, 334 , 330 , 351 , 347 , 368 , 364 , 385 , 381 , 402 , 398 
419,415, 436 , 432 , 453 , 449 , 470 , 466 , 487 , 483 , 504 , 500 
1 7 , 34 , 51 , 68 , 85 , 1 02 , 1 1 9 , 1 36 , 1 53 , 1 70 , 1 87 , 204 , 221 , 238 , 255 , 272 
289 , 306 , 323 , 340 , 357 , 374 ,391,408, 425 , 442 , 459 , 476 , 493 

96. 1  , 17 

198.1.17 

300.1.17 

402.1.17 

■End  o-f  Problem" 

3,0.0001 


# 
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769 


SAMPLE  OUTPUT 


CE504CW 


Userid:  CE504CW 

File:  CE504CW  OUTPUT 


USERID 

CE504CW 

ORIGIN 

RSCS 

DISTCODE 

CE504CW 

SYSTEM 

GITVM1 

SPOOLID 

8001 

RECORDS 

00001063 

FILE 

CE504CW 

OUTPUT 

CREATED 

07/03/89 

23:02:08 

PRINTED 

07/03/89 

23:02:23 

CLASS 

A 

FORMS 

XD 

DEVICE 

600 

SEQUENCE 

769 

&?SS2\ 


F  ,S  lot. 

1  .  IOS 

-L  .  0 

2  ,0101 

<4  .  OOZG 
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F  Kit  It 


se«o  i 


WIFE  .  16 - HIIIHC  .  3J 

NOOK 

NOOK 

NOOK 

NODE 

1 

3 

4 

7 

uatnr 

* 

NOOK 

1  1 

• 

NOOK 

It 

NOOK 

20 

Udot 

'  a . 

NOOK 

24 

NOOK 

24 

NODE 

24 

"•'  *  46  6006 
-i»o. 0000 
-1*0.  0000 

-  HO  OOOO 


•O . OOOO 
1 « . OOOO 
14 . OOOO 

T4"*666o 

I  4  . OOOO 
I  4  . OOOO 

I  4  OOOO 


— ¥V*r6'6" 

1  1  *431 

1*  .  3  1  47 
32  4*00 


34 . 4400 
3 . 4400 
I  .  44  J4 

1  4  4*00 

2* . 4400 
32  4*00 


4  4  4  Z  4 

12  .  1*42 

13  4  2*4 

ilKlH 

24 . 4244 
.  4440 
4.4443 


1 4 . 4440 

22 . 1 440 

27 . 4440 

"oVobob" 
4  .  3333 
10.4*47 
1 4  OOOO 


# 

NOOK 

NOOK 

NOOK 

7* 

4  4 

44 

mot 

46 

NOOK 

42 

NOOK 

*4 

NOOK 

*4 

0  .  OOOO 
IO . OOOO 

1 0 . 0000 


1  .  4000 
1  .  4000 
1  .  4000 

Tl'iobb' 

1  .  4000 
4  .  4  000 
«  4000 


[.a  f.T.T.i 


24 . OOOO 
•  .  3400 
4  1 0OO 

"ib .  4466" 

1  4  .  OOOO 
1  4  .  OOOO 
20 . OOOO 


4 . 0440 
10  4321 
1 4  OOOO 

"ft  Vooob" 

20 . OOOO 
• . 4774 
4 . 4443 


24 . OOOO 
24  OOOO 
24 .OOOO 

2  0>.  OOOO 

3  1  .  4000 
31 . 4000 


NOOK  177 
HOOK  174 
NOOK  141 

•  MttttfUi' . 

NOOK  1*0 
NOOK  1*3 
NOOK  1*4 


NOOK  1*4 
NOOK  204 
NOOK  207 

mM'M . 

NOOK  211 
NOOK  213 
NODR  214 


111»H 


OOK  224 
000  22* 
OOK  324 

OOK  232 
OOK  23* 
OOK  241 


NOOK  344 
NOOK  347 
NOOK  344 

iiOttt'"*** . 

NOOK  244 
NOOK  240 
NOOK  342 


NOOK  244 
NOOK  373 
NOOK  274 
■****"•• 
NOOK  27* 
NOOK  2*1 
NOOK  243 


44 . 400 
4* . 400 

44 . 400 

*t:  too 

*1 . 400 

4 1 . 400 
*  1  . 400 


1  * . OOOO 
1  *  . OOOO 
22 . OOOO 
"‘•"■*64©" 

4.7447 
10 . 34*3 

1 4 . OOOO 


10 . 37*0 
14 . OOOO 
14 .OOOO 
22.0000 
•  1  .  1024 
4.4*43 


14.0000 
f  * . OOOO 
22  .  OOOO 

•■•i'.’imr 

4.44  17 
IO . 2404 
1  *  .  OOOO 


20 . OOOO 
-  1 . 3300 
4 . 44*7 

TOViim  ■■ 

1  * .  0000 
14 . 0000 
20 . 0000 


4 . 3300 

10 . 1 *40 

1  *  .  OOOO 

"WVbbttb"' 

22 . OOOO 
•  1 . 732* 

4 . 17*3 


1  I . OOOO 
1* . OOOO 
22 . OOOO 

VV*624" 

4 . 1317 
10.04*4 
1  * . OOOO 


NOOK 

292  K 

• 

NOOK 

2*4  1 

NOOK 

2*4  X 

, 

Ndttt 

NOOK 

300  I 

NOOK 

307  X 

NOOK 

30*  I 

NOOK 

313  I 

NOOK 

314  1 

NOOK 

317  K 

m . * 

NOOK 

32*  X 

• 

NOOK 

334  X 

NOOK 

330  X 

NOOK  334 
HOOK  34 1 


4 . 4000 
4 . 4000 
4 . 4000 
4  4666 

*  4000 

*  4000 
f  4000 


22 . OOOO 
•  2 . 0300 
3 . *400 

"’4 ’.■*■*  66"' 

1  *  OOOO 
1  * . OOOO 
22 . OOOO 


3 . 4433 

• .*317 
1  * . OOOO 

T47 6666 

20 . OOOO 
-  2  247* 
3 . *2*3 


I  4 . OOOO 
1 4 . OOOO 
20  OOOO 

■■•2'V*a'24"' 

3  7  117 

4  .  •**• 

1  *  .  OOOO 


•I  222  233  231  240  241 


22  1  2  1  *  2  14 

234  224  223 

??«  27* 


200001*00 

20000* *00 


. 290001*00 
. 200001*00 

20000*  *00 


2*7 

2  ft  • 

2  ft  ft  2  ft  ft 

2  ft  ft 

2  ft  ft 

24ft 

>«• 

. 200001*00 

2S  ft 

2  ft  7 

273  77* 

27ft 

2  ft  ft 

2  ft  ft 

2  ft  7 

.  200001*00 

2  ft  ft 

2  ft  ft 

2  7  ft  27ft 

277 

2  ft  ft 

2  to 

2  ft  ft 

. 200001*00 

.  200001*00 

2  fto 

2«ft 

277  27ft 

27  ft 

27o 

212 

2ft  1 

. 200001 *00 

. 266601*66 

ESI 

Kfil 

2ft  1 

27  1 

2  ft  ft 

2  ft  3 

. 20000C *O0 

(  VWiiHII  IWI 1 


I  1 1 


2  f  ft  2  ft  ft 

2  7  ft  27ft 


.  2 0000ft *00 
2ooooe*oo 


.  200001 *00 
.  20O00C  «00 


•  a 

277 

a«a 

2*4 

2ft  ft 

2ft  ft 

247 

27* 

27ft 

20000C  *00 

.200001*00 

•  4 

27ft 

2  ft  7 

2ft  ft 

2ft  7 

2  V  ft 

2 1  ft 

2ftO 

. 200001*00 

.  20O00C *00 

• 

it 

2ft  1 

2  ftft 

2  ft  ft 

2  •  • 

300 

2  ft  • 

2  ft  3 

2  ft  2 

. 200001*00 

.  20000C*00 

ft  ft 

•  ‘}ft  o 

Sol 

567  " 

•  564  " 

“36  ft" 

362" 

"24  4" 

"26'f . 

. .2 66664*66 . 

.266606*66 

•  7 

2ft  2 

302 

sot 

3  10 

a  1 1 

303 

2  ft  4 

2  ft  3 

. 260001*00 

.200001*00 

•  ft 

2*4 

303 

3  1  1 

3  1  2 

313 

304 

sac 

2ft  ft 

. 200001*00 

.200001*00 

•  ft 

2  ft  ■ 

304 

3  1  3 

3  1  4 

3  1  ft 

soft 

3  •  ft 

2ft  7 

. 200001*00 

200001*00 

ft  1 
•  2 
•  3 

307 

304 

3  1  1 

31ft 

31ft 

320 

324 

32ft 

32ft 

32S 

327 

32ft 

32ft 

32ft 

330 

31ft 
320 
32  1 

30ft 

3  1  1 

a  i  a 

30ft 

3  10 

312 

. 200001*00 
. 200001*00 
. 200001*00 

. 200001*00 

ft* 

4  » 3 

22  i 

35© 

33) 

332 

343" 

"i  14" 

’  3l4 . 

. 26o6ot  *66 

. 366661*66 

•  ft 

31ft 

322 

332 

333 

334 

323 

3  1  7 

3  1  ft 

. 200001*00 

•  ft 

32« 

33S 

34  1 

342 

343 

33ft 

33ft 

32ft 

. 200001*00 

.  200061*00 

•  7 

Em 

33ft 

343 

344 

34  ft 

337 

321 

327 

. 200001*00 

. 200001*00 

347 

34ft 

34ft 

334 

333 

33  1 

. 200001*00 

. 200006*00 

33ft 

34  i 

3ft© 

3S  1 

340 

334 

333 

. 200006*00 

. 200006*00 

3*2 

314 

Sftft 

3.0 

343 

343 

342 

. 200006*00 

. 200001*00 

’343 

“366“ 

"Sit 

"342 

“364" 

"346" 

"544 . 

.266666*66 . 

""'■'266666*60 . 

3  ft  4 

3  C  2 

J  ft  3 

3*4 

311 

347 

34ft 

. 200006*00 

. 200006*00 

SIS 

3  44 

ass 

Sftft 

a  s  4 

34ft 

34ft 

. 200006*00 

. 200006*00 

3ft  ft 

3  ftft 

3  ft  7 

>ni3 

3  5  7 

39  1 

3  ft  0 

200006*00 

. 200006*00 

107 

3  ft  0 

370 

377 

37ft 

37ft 

371 

3  «2 

3ft  1 

. 200006*00 

. 200006*00 

ioi 

a  ft  2 

371 

37ft 

3  ftO 

3ft  1 

372 

3  8  4 

3*3 

. 200006*00 

.  200006*00 

1  Oft 

3«4 

322 

3ft  1 

3  ft  2 

3*3 

373 

Sftft 

3  ft  ft 

200006*00 

.  200006*00 

# 

1  1  o 

"364’" 

373 

563 

3  ft  4 

374 

"3  64“ 

"567 . 

. 2 66606 *66 
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Using  a  Personal  Computer  in  Running 
the  Aral  Seepage  Program 

I .  Introduction: 

The  Aral  Seepage  Program  is  located  on  the  Georgia 
Institute  of  Technology  main  frame  computer  in  the  Cyber  B 
area.  To  run  many  runs  of  the  program  it  is  very 
convenient  to  prepare  and  edit  the  data  files  on  a 
personal  computer,  using  an  editor  or  word  processor.  A 
brief  description  of  the  method  for  this  particular 
program  follows.  A  user  should  get  instruction  and 
information  for  use  of  the  Cyber  B  from  the  Office  of 
Computer  Services,  Rich  Building,  Georgia  Institute  of 
Technology. 

II .  Preparing  the  Data  Files: 

Using  the  Aral  Seepage  Program  User's  Manual  (Pirtle, 
Appendix  A)  draw  and  label  the  mesh  for  the  problem  to  be 
run.  An  example  is  Figure  3-5  of  this  study,  but  a 
working  copy  of  the  mesh  should  be  large  enough  for  easy 
referral  and  use.  With  the  mesh  and  a  complete  list  of 
problem  characteristics  a  data  file  can  be  prepared,  from 
the  user's  manual  instructions. 

In  this  study  the  number  of  nodes  and  basic  shape  of 
the  mesh  was  to  be  unchanged  for  all  of  the  runs.  To 
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simplify  preparing  the  data  files  a  spreadsheet  was 
prepared.  Changing  the  dimensions  or  hydraulic 
conductivities  (vertical  or  horizontal)  required  one  entry 
and  the  spreadsheet  adjusted  the  node  positions  from 
simple  formulas.  The  complexity  of  the  spreadsheet  is 
entirely  dependent  on  the  problem  considered  and  the  skill 
of  the  user  with  a  spreadsheet  program.  Any  of  the 
commercially  available  programs  will  work  satisfactorily. 
If  the  study  will  involve  several  different  mesh 
arrangements  and/or  only  a  few  runs  the  data  file  can  more 
easily  be  prepared  on  an  editor  by  typing  it  in,  if  the 
editor  is  set  to  create  an  ASCII  file. 

An  ASCII  file  is  a  file  of  the  characters  as  entered. 
Word  processing  programs  and  spreadsheet  programs  have 
hidden  command  codes  which  control  the  programs  so  that 
the  features  are  executed.  Features  are  the  calculation 
of  values  in  a  spreadsheet  or  the  margins,  justification 
and  format  of  a  word  processor. 

Most  spreadsheet  programs  have  a  function  or 
ancillary  utility  program  for  conversion  of  the 
spreadsheet  to  a  values  separated  by  commas  ASCII  format. 
With  that  utility  the  spreadsheet  no  longer  will  do 
calculations  or  have  the  simple  column  alignment. 

In  this  study  the  programs  available  provided  an 
ASCII  file  with  all  values  separated  by  commas.  The  file 
did  have  commas  to  separate  what  had  been  empty  cells 
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within  the  spreadsheet.  Other  programs  may  not  have  that 
result.  To  delete  the  commas  a  word  processor  was  used. 
A  program  which  allows  editing  an  ASCII  file  is  necessary 
for  this  operation.  If  a  word  processor  is  used  in  the 
normal  mode  of  edit  command  codes  will  be  imbedded  to  set 
margins  and  the  file  will  not  run  properly  on  the  Aral 
Seepage  Program. 

Ill .  Sending  the  Data  File: 

With  a  proper  data  file  prepared  running  the  Aral 
Seepage  Program  from  a  personal  computer  requires  one 
device  and  several  clearances.  An  authorization  to  use 
the  Cyber  B  must  be  obtained,  with  a  file  storage  area, 
a  password  and  an  identification  code.  This  is  explained 
by  manuals  available  at  the  Office  of  Computer  Services. 
Access  to  the  Aral  Seepage  Program,  named  OWRCC  must  now 
be  obtained  from  the  School  of  Civil  Engineering  at  the 
Mason  Building.  The  device  needed  is  a  MODEM  for  transfer 
of  data  and  commands  to  the  Cyber  B.  The  speed  and 
manufacture  of  the  MODEM  is  subject  to  the  needs  of  the 
user. 

A  communications  program  is  necessary  to  use  the 
MODEM  and  many  such  programs  are  available.  Georgia  Tech 
offers  students  and  faculty  a  program  named  COGITATE  for 
use  in  communicating.  In  this  study  a  different  program 
was  used  but,  it  had  the  KERMTT  file  transfer  program  and 
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this  is  necessary. 

With  the  communications  program  call  the  GTNET  system 
and  when  a  connection  is  made  enter  the  GTNET  command: 

>>  Connect  *OCS  CYBER  B 

A  connection  to  this  directory  will  be  made  and  right  to 
access  will  be  challenged: 

Family  CYBERB 

Userid  (input  the  assigned  ID) 

Password  (input  the  selected  password) 

Access  will  be  granted  to  the  Cyber  B  area.  To  copy  the 
data  file  to  your  file  for  execution  of  the  Aral  Seepage 
Program  the  program  COGITATE  has  menu  entries  to  execute. 
If  another  communication  program  is  used  the  KERMIT 
protocol  is  necessary.  The  details  of  this  description 
are  better  left  to  someone  else.  In  this  study  the 
following  procedure  worked. 

\  KERMIT2 

Control  was  now  with  the  user’s  computer  by  use  of  KERMIT 
commands.  Executing  the  ’’Send"  command  of  KERMIT  and 
naming  the  data  file  to  be  downloaded  successfully- 
transferred  the  file.  The  file  name  was  changed  as  the 
cyber  will  use  the  first  eight  alpha-numeric  characters 
of  the  name  and  ignore  any  periods  normally  used  to 
separate  suffixes  in  the  MS-DOS  environment.  At  the 
completion  of  sending  data  files  execute  the  KERMIT 
command  to  "Finish"  and  control  will  be  re-established 
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with  the  Cyber. 


IV.  Running  the  Program: 

To  run  the  program  and  receive  a  printed  output  from 
the  laser  printers  in  the  Rich  Building  do  the  following: 


\  OLD,OWRCC 
\  GET, DATA 


(note: DATA  is  data  file  as 


named  upon  receipt  by  Cyber.  To  see  names 


answer  \  with  LF) 

\  OWRCC , DATA , PRINT 


(note: PRINT  is  name 


selected  for  the  output  by  user.  If  no  print 
is  specified  the  output  of  the  run  will  go  to 
the  screen.  With  a  MODEM  this  is  a  very  long 
process  and  is  not  recommended.) 

A  short  wait  of  less  then  a  minute  can  normally  be 
expected  here.  If  an  error  message  is  received  the  data 
file  should  be  checked  for  errors  in  concept  and  tried 
again.  If  a  transfer  error  is  suspected  send  the  file  a 


second  time. 


\  SAVE, PRINT 


(note: SAVE  only  if  a  copy  of 


the  output  is  desired  for  the  holding  on  Cyber) 
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\  PRINT?  (note: A  hard  copy  of  the 

program  output  will  be  received  in  the  Rich 
Building  if  the  output  file  name  PRINT  is  input 
at  the  first  query  and  all  other  entries  are 
allowed  to  default  by  responding  with  the 
“ENTER"  or  "RETURN"  key.  The  last  query  asks 
for  the  USERID  and  responding  with  the  ID  has 
it  printed  on  the  cover  of  the  hard  copy  for 
easier  identification. 

There  are  several  other  printer  locations  and  many 
options  one  can  exercise  in  printing  the  results.  For 
more  information  consult  the  advisors  in  the  Rich 
Building. 

To  exit  from  the  Cyber: 

\  LOGOUT 
»  LOGOUT 
V.  Caution: 

This  appendix  is  not  a  replacement  for  the  user's 
manuals  of  the  programs  used,  nor  for  the  instructions 
and  guidance  of  the  Office  of  Computer  Services.  It  is 
a  quick  and  easy  guide  of  a  procedure  that  has  been 
successful  in  running  programs  on  the  Aral  Seepage 
Program. 


APPENDIX  C 


6 


# 

I 

As  stressed  above  many  different  programs  will  do 
%  the  functions  of  the  above  procedure.  For  information 

the  following  were  used: 

Spreadsheet:  SuperCalc  3 

•  Editor:  WordStar  4.0 

Communications:  ProComm  2.1 

It  would  be  an  advantage  if  the  output  file  from  the 

•  Aral  Seepage  Program  could  be  captured  and  downloaded 
from  the  Cyber.  This  was  attempted  with  KERMIT  once  but, 
was  discontinued  after  twenty  minutes  of  transfer  with  a 

•  1200  Baud  MODEM. 


# 
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Arbitrary  Selection  of  Radius 

In  this  study  the  radius  of  influence  was  determined 
as  the  distance  at  which  no  more  then  a  10%  decrease  in 
total  flow  over  a  100  foot  change  in  trial  radii  occurs. 
This  arbitrary  selection  is  certainly  open  to  discussion. 
It  is  very  common  for  less  restrictive  criteria  to  be  used 
in  the  definition  of  radius  of  influence  of  drawdown. 

Figures  E-l  and  E-2  are  cross  section  drawings  of 
the  free  water  surface  in  the  trench  area  if  the  radius 
of  influence  had  been  selected  at  lower  values  based  on 
some  different  arbitrary  standard.  Obviously  the  free 
surface  reaches  the  upper  blanket  and  possibly  the 
foundation.  The  formulas  used  in  Chapter  4  were  based  on 
field  data  of  seepage  and  the  field  interpretation  of  the 
radius  of  influence,  under  some  standard  determination. 

If  the  selection  of  radius  of  influence  for  this 
study  had  been  based  on  a  different  criteria  the 
predictive  forms  and  dimensionless  products  would  have 
been  changed.  Using  a  standard  of  0.1  ft5/day/ft,  maximum 
change  in  trial  radii  of  influence,  Criteria  A,  results 
in  Figures  E-3  and  E-4.  These  figures  are  the 
dimensionless  product  for  a  different  standard  of  radius 
of  influence.  The  new  standard  is  on  average  0.07 
ft3/day/ft  less  restrictive  and  predict  the  radii  of 
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Table  E-l.  The  new  curves  of  dimensionless  products  in 
Figures  E-3  and  E-4  show  the  less  restrictive  curves  to 
be  nearly  unchanged. 

As  a  further  check  a  second  standard,  Criteria  B  was 
developed.  This  standard  is  defined  as  a  total  flow 
increase  of  75%  from  the  total  flow  at  a  trial  radius  of 
700  feet.  The  scatter  plot  of  the  dimensionless  product 
versus  total  flow  and  maximum  free  surface  is  shown  by 
Figures  E-5  and  E-6  for  Criteria  A  and  Criteria  B.  Table 
E-7  summarizes  the  results  of  Criteria  B. 

Figures  E-7  and  E-8  show  the  final  curves  for  total 
flow  and  the  maximum  free  surface  against  the  dimension¬ 
less  product  for  all  considered  standards.  The  curves 
are  little  changed.  The  use  of  one  composit  curve  would 
suffice  in  predicting  total  flow  and  maximum  free  surface 
for  any  standard  of  radii  determination. 

An  incorrect  radius  of  influence  will  under  predict 
the  flow  and  maximum  free  surface  if  too  high.  A  con¬ 
servative  approach  would  be  to  reduce  the  radius  of 
influence  as  predicted  from  Figures  4-1,2  and  3.  The 
amount  of  adjustment  would  require  additional  field  data 
from  a  trench  drain  system. 
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FREE  SURFACE  PROFILE 

Slope  =  0%,  L  =  200',  Spacing  = 

H  =  60',  h  =  20’ 
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Free  surface  between  trenches 


TABLE  E-2 


SUMMARY  OF  DATA,  CRITERIA  B  RADIUS  DETERMINATION 
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TOTAL  FLOW 


Radius  selected  from  alternate  method 


Selection  of  L  from  alternate  method 
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Fig  E-5 

Comparison  of  two  alternate 
criteria  for  radii  determination 
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Total  Flow 
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Fig.  E-7 

Primary  and  two  alternate 
radii  determination  methods 


Maximum  Free  Surface 

Dimensionless  Comparison 
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Fig.  E-8 

Primary  and  two  alternate 
radii  determination  methods 


